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INTRODUCTION 

Makers of atmospheric models could well 
be intimidated by the wide variety and unusu- 
al character of the red-giant stars as revealed 
by the observations already described, which 
quite commonly include such features as emis- 
sion lines, thermal emission from dust shells, 
maser emission, mass outflow (and inflow), po- 
larization, variability, and rapid evolution. 
These phenomena present tough questions. To 
what extent are the atmospheres of M, S, and 
C stars to be understood simply as cooler ex- 
amples of ordinary warmer stars? Are inhomo- 
geneous atmospheres common among red gi- 
ants and supergiants? Will discrepancies due to 
our present inadequate treatments of convec- 
tion conceal other theoretical deficiencies? Are 
chromospheres permanent or transient? By 
what mechanisms are they produced? Where 
and how does grain formation and growth oc- 
cur? What mechanisms are responsible for mass 
loss, and what are the consequences? With ex- 
tended atmospheres and temperature inver- 
sions, aren't departures from local thermody- 
namic equilibrium (LTE) likely? Can one neatly 
separate atmospheres into photospheres, chro- 
mospheres, and circumstellar shells? The ques- 
tions go on and on. 

In spite of the uncertainties raised by these 
questions, atmospheric model makers have his- 
torically treated all these more exotic phe- 


nomena as secondary and have concentrated 
first on modeling photospheres. A photospheric 
model is not only a desirable first step, but it 
has a rather logical priority since all other struc- 
tures, processes, and phenomena depend on it. 
In this chapter, we describe the static thermal 
atmosphere and compare it$ predictions to ob- 
servations both to test the validity of the classic 
assumptions and to distinguish and describe 
those spectral features with diagnostic value. 
By thermal atmosphere, we mean an atmo- 
sphere based on hydrostatic equilibrium, radia- 
tive equilibrium, and LTE. Nonthermal atmo- 
spheric models, including nonradiative heating 
and departures from LTE, are treated in a suc- 
ceeding chapter (de la Reza, this volume), and 
more exotic phenomena are described in other 
chapters. 

Several useful reviews of the challenges of 
computing atmospheres for red-giant stars and 
using these to interpret observations have al- 
ready been given. Vardya (1970) gave a broad- 
brush treatment of the topic from information 
then available, and Johnson (1972) analyzed the 
few theoretical models computed to that time. 
Although advances in modeling atmospheres of 
red-giant stars have come slowly, it is hearten- 
ing to note the obvious progress from these 
early reviews to the present. Remarkably, even 
the recent excellent general works on stellar at- 
mospheres (for example, Gray, 1976; Mihalas, 
1978; Baschek and Scholz, 1982) contain few 


323 



applications to red-giant stars. More valuable 
are recent specialized reviews. Carbon (1979) 
gives a detailed genera) Review of problems and 
techniques of computation of atmospheres of 
intermediate and late-type stars from the point 
of view of the model maker; Gustafsson (1981) 
writes a topical review of late-type stars from 
the point of view of the observer; and Johnson 
(1985) discusses the current availability of mod- 
els for peculiar red-giant stars. Carbon (1984) 
examines both the opacity-distribution function 
and the opacity-sampling treatments of line 
opacities, the latter by means of a Monte Carlo 
technique. 

PHYSICAL PRINCIPLES 

For modeling red-giant atmospheres, work- 
ers have adopted the same principles and tech- 
niques used in modeling hotter main-sequence 
stars, for which these principles were perhaps 
more appropriate. The most obvious exception 
is in the treatment of line opacity, for which 
the relatively few lines in hot stars can be 
treated individually but the millions of lines in 
cooler stars must be treated collectively. The 
basic physical principles are hydrostatic 
equilibrium, radiative equilibrium, and local 
thermodynamic equilibrium (LTE). The solu- 
tion of these equations in a plane-parallel 
horizontally homogeneous geometry constitutes 
the problem of the classic stellar atmosphere. 

Hydrostatic equilibrium in a plane-parallel 
geometry is: 

dp = -ggdz , (7-1) 

where the total pressure p = p(gas) + p(radia- 
tion) + p(turbulence), and the surface gravity 
g = GM/r 2 is held constant. Turbulent pres- 
sure is often neglected. 

The constancy of the total energy flux is: 

F = a r* ff , (7-2) 

where F is the energy flux (erg cm 2 s' 1 ) and 
r is the effective temperature, defined by 

eff . 

Equation (7-2). Strict radiative equilibrium 1 m- 


plies j F v dv = F(rad) = F, whereas convec- 
tive energy transport is included through F(rad) 
+ F(conv) = F. When included, convection 
has usually been treated by the usual local mix- 
ing length (LML) theory without overshoot 
(Bohm-Vitense, 1958; Henyey et al., 1965; Mi- 
halas, 1978; Lester et al., 1982). 

As generally used, LTE includes two key as- 
sumptions: (1) concentrations of molecules, 
atoms, and ions are given by the relevant equa- 
tions of equilibrium statistical mechanics with 
the single parameter of electron temperature; 
and (2) the source function for all bound-bound 
transitions and for all continuous transitions ex- 
cept for pure scattering is the blackbody radia- 
tion (Planck) function, B{v, T). Nothing in the 
modern treatments of bound-bound opacity 
(opacity distribution functions, opacity samp- 
ling (OS), or Voigt-analog/Elsassar-band 
model; see following section) requires this 
assumption on the line source function; they are 
all flexible enough to accommodate a simple 
mixture of pure absorption and pure scatter- 
ing. (The OS has, however, a clear advantage.) 
The LTE assumption is made for simplicity and 
to test this hypothesis in the face of departures 
anticipated because of the observed departures 
of the emitted energy flux of cool stars from 
that of a blackbody and the low gas densities 
(hence collision rates) in the outer atmosphere. 

We attempt to display the present situation 
and foreseeable developments in Figure 7-1. 
Principles of the classic stellar atmosphere are 
shown at the left. Developments or improve- 
ments currently available or on the horizon are 
shown in solid boxes; developments still in the 
future are shown in dotted boxes. Except for 
the addition of convective energy transport and 
the extension to spherical geometry, both of 
which are, in some form, already in use, gen- 
eralizations of these classic principles may be 
slow in appearing. Recall that we are speaking 
here of stellar photospheres. Some of the im- 
provements labeled “future” are regularly ap- 
plied in studies of circumstellar (CS) envelopes 
and mass loss; current progress in these areas 


324 



HE = hydrostatic equilibrium 

RE = radiative equilibrium 

LTE = local thermodynamic equilibrium 

PPG - plane-parallel geometry 

SG = spherical geometry 

CE = convective equilibrium 

SSE = steady-state equilibrium 


Figure 7-1. Overview of the physical principles 
underlying current theoretical research on clas- 
sic thermal atmospheres of red-giant stars. 
Present physical principles and generalizations 
are enclosed with solid lines , and possible future 
generalizations are enclosed with dotted lines. 


can be gaged from other chapters in this vol- 
ume. The rate of future progress can be guessed 
at by the reader as it has been guessed at by 
others (Johnson, 1985). With the arrival of in- 
formation in the ultraviolet, infrared, and 
radio regions and with the prospects for em- 
ploying supercomputers to model the entire at- 
mosphere, including mass flow, on the horizon, 
we appear to be entering an exciting era. 

OPACITIES 
Continuous Opacities 

Many of the opacities for cool stars are 
similar to, or can be extrapolated from, those 


at higher temperatures. Throughout most of the 
photosphere, H" provides the dominant con- 
tinuous opacity over most of the spectrum even 
down to stars of T ff = 2500 K and for a 
variety of chemical compositions (Johnson, 
1982). In the coolest region of the atmosphere, 
hydrogen may be associated into H 2 , and the 
H“ opacity will become so small that He" (free- 
free) will become a dominant continuous opac- 
ity. This may be the case in hydrogen-deficient 
atmospheres as well. Bound-free continua of 
H I are only of minor importance, and the 
bound-free continua of He I are completely 
negligible. (These may be of importance, how- 
ever, in shocks or in the hot chromospheres of 
cool stars.) Absorption due to bound-free con- 
tinua of other neutral elements (Si, Al, Na, Mg, 
Ca, and Fe) are likewise of minor importance 
because their absorption edges lie in the ultra- 
violet, where the stellar flux is weak. Except for 
H", the opacity of negative ions (Vardya, 
1970) does not appear to be significant. Al- 
though Rayleigh scattering from H I, H 2 , and 
He I — with minor contributions from C I, N 
I, O I, and N 2 — is most significant in the 
ultraviolet, its effect extends increasingly long- 
ward in wavelength as the temperature de- 
creases, and it becomes of major importance 
in the upper photosphere. Information on cross 
sections for continuous opacity sources is con- 
tained in, for example, Carbon and Gingerich 
(1969); Kurucz (1970); Tsuji (1971); and Bas- 
chek and Scholz (1982); while references may 
be found in Gustafsson et al. (1975); Mihalas 
(1978); Kurucz (1979); Johnson et al. (1980); 
and Dzervitis (1983). 


Dust Grains 


Wherever condensation occurs in a stellar 
atmosphere, dust will almost certainly become 
a major opacity because of its very large ab- 
sorption coefficient. Many cool stars, especially 
cool supergiants and Miras, show evidence of 
circumstellar dust (Lefevre, this volume). 
Whether grains commonly form in the outer 
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stellar photosphere is still uncertain, despite 
considerable work on this important problem. 
Two enormous obstacles to progress in over- 
coming this uncertainty are the lack of knowl- 
edge of both the physical conditions in the outer 
photospheres (this chapter; Zuckerman; 1980; 
Goldberg, this volume) and the mechanism of 
nucleation (Deguchi, 1980; Draine, 1981; 
McCabe, 1982; Alexander et al., 1983; Gail and 
Sedlmayr, 1984). Vital to the question of grains 
in red-giant atmospheres is the location of the 
condensation point, but this is still quite uncer- 
tain (c.f. Draine, 1984). For clean silicate grains 
in M stars, condensation may occur at a dis- 
tance (from the stellar center) of 3 to 4 stellar 
radii (Deguchi, 1980) or perhaps closer (Draine, 
1981; Alexander et al., 1983). In carbon-rich 
atmospheres, condensation of graphite may oc- 
cur at 5 stellar radii (from the stellar center ), 
and SiC may condense even closer (McCabe, 
1982), perhaps as close as 1 .5 stellar radii (Lucy, 
1976; Woodrow and Auman, 1982). Several 
minerals condense at about 1500 K (cf. Draine, 
1981; McCabe, 1982; Alexander et al., 1983), 
and such temperatures might be reached in the 
outer photospheres of sufficiently cool stars, 
especially if presently unaccounted for spherical 
effects (Schmid-Burgk and Scholz, 1981) lower 
the temperature even further than predicted by 
current models. On the other hand, chromo- 
spheric heating may raise the temperatures in 
the outer layers. Clearly, the most favorable site 
for condensation is the outer atmosphere of a 
cool Mira variable star (Deguchi, 1980; Alex- 
ander et al., 1983), where grain formation may 
be related not only to mass loss but also to 
pulsation (Woodrow and Auman, 1982). At 
present, there is no compelling evidence for dust 
in any cool-star photosphere (see, for example, 
Chapters 3 and 7 and the foregoing references). 
For this reason, and for the sake of simplicity, 
dust is neglected in most model atmospheres. 
If future observations demonstrate the existence 
of photospheric dust grains, their inclusion will 
extend atmospheric modeling into a rather dif- 
ferent opacity regime. 


Spectral Line Blanketing 

Before plunging into the complexities of the 
subject, it is worth reminding ourselves of the 
central role of the atomic and molecular data 
themselves. Much work lies ahead, particularly 
for the molecules, and we salute those who pur- 
sue this unglamorous but necessary research. 

Several excellent reviews of line blanketing 
are available. A description of LTE line blan- 
keting in terms of physical mechanisms is given 
by Carbon (1979), and a mathematical analysis 
is presented by Mihalas (1978). Quantitative 
non-LTE analyses are provided by Athay and 
Skumanich (1969), Mihalas and Luebke (1971), 
and Athay (1972). Much of this development 
builds on foundations laid earlier by Chandra- 
sekhar (1935), Munch (1946), and Pecker 
(1951). 

From the point of view of the observer, the 
most striking feature of a stellar spectrum is the 
removal of light by the spectral lines, an effect 
referred to as line blocking . Since the line 
represents an increase of opacity and the total 
flux is required to be constant in radiative equi- 
librium, the temperature in deeper layers must 
rise to drive the same flux through the increased 
opacity, an effect referred to as backwarming. 
In addition, the gas in the outer layers may be 
heated or, more frequently, cooled. Altogether, 
the effects are referred to as line blanketing. 

Physical processes near the surface and their 
effect on the thermal structure are complex. 
Purely coherent scattering lines change the tem- 
perature very little because they are weakly 
coupled to the thermal reservoir of the atmo- 
sphere (Chandrasekhar, 1935; Munch, 1946; 
Mihalas, 1978). Purely absorbing lines generally 
act to lower the temperature of the outer layers 
if they lie on the redward side of the flux max- 
imum (e.g., H 2 0, CN, and CO); purely ab- 
sorbing lines lying blueward of the flux peak 
(TiO and most atomic lines) tend to raise the 
temperature of the layers in which they finally 
become thin. These latter effects can be under- 
stood as follows (Krupp et al., 1978; Carbon, 
1979; Gustafsson, 1981). Absorption of a red 
(low-energy) photon by an atom or molecule 
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and the reemission of a (higher energy) photon 
characteristic of the local temperature in the 
outer layers results in a net energy loss (i.e., a 
cooling). In the same way, absorption of a UV 
photon and reemission of a (lower energy) pho- 
ton characteristic of the local temperature re- 
sults in a net energy gain (i.e., a heating). This 
is a simplified picture of a very complex pro- 
cess, however, and a careful quantitative de- 
scription, including the effects of departures 
from LTE and a proper description of the 
mechanism of line formation may be necessary 
to understand the physics involved in a partic- 
ular line (Jefferies, 1968; Athay and Skuma- 
nich, 1969; Mihalas and Luebke, 1971; Athay, 
1972; Mihalas, 1978). Although line formation 
in the general case, where LTE does not hold, 
properly belongs in a chapter on nonthermal 
phenomena (de la Reza, this volume), we men- 
tion it here because of its possible importance 
in treating line blanketing. 

In fact, the effect of a line on the atmo- 
spheric structure depends on the line strength 
relative to the continuum, the variation of the 
line strength with depth, the wavelength of the 
line, and the mechanism of line formation. This 
last consideration is sometimes specified before- 
hand (e.g., pure absorption or pure scattering), 
which oversimplifies the actual physics of the 
problem. Of these effects, the blocking and 
backwarming depend principally on the line 
strength, whereas the surface cooling depends 
principally on the mechanism of line formation, 
but also on the line strength (Mihalas and 
Luebke, 1971). Quite generally (cf. Athay and 
Skumanich, 1969; Athay, 1972), the condition 
of radiative equilibrium requires dH/dr = 0, 
where H is the normalized total (integrated) 
flux. Although the total flux is conserved, the 
continuous flux is now no longer conserved but 
is given by: 


<W C M - 1 

dr r 

C O 


+ 



-J v )dv 


t + c , 


(7-3) 


where M = Av D (for a Doppler line); Av D 
is the Doppler halfwidth; Z = 1 - J/Sj is the 
net radiative bracket; J is the mean intensity 
averaged over the absorption coefficient; and 




(7-4) 


Here £ and c refer to line and continuous quan- 
tities, respectively, and v designates total (line 
plus continuous) quantities. The net radiative 
bracket (Thomas, 1960) or escape coefficient 
(Athay, 1981), when multiplied by the Einstein 
spontaneous emission coefficient, gives the net 
downward rate of radiative transitions in a 
given spectral line. This equation states that the 
presence of a spectral line causes a net exchange 
of energy between the line and the continuum— 
an exchange which may proceed in either direc- 
tion and may depend sensitively on depth. It 
is the atmospheric adjustment to the re- 
quirements that total flux be conserved which 
produces the changes in the thermal structure 
referred to as line blanketing. 

The c term (Athay, 1972) is proportional to 
the equivalent width of the line and is defined 
to be positive if the line is in absorption. This 
term is then independent of the mechanism of 
line formation, at least to the first order. If an 
absorption line is present, it provides a sink for 
the continuous photons and thus tends to cool 
the atmosphere. This is always the case in the 
outer layers, and c therefore acts to cool the 
outer layers regardless of its position in the 
spectrum, strength, or depth of formation. As 
one proceeds into the atmosphere, the line may 
stay in absorption or may go into emission, de- 
pending on the wavelength, the line strength (or 
value of r o ), and the depth dependence of r . 
Lines in the blue part of the spectrum— blue- 
ward of the flux peak, where J > B — tend to 
stay in absorption and therefore cool (provide 
a sink for continuum photons) at all depths, 
whereas lines redward of the peak, where J < 
B, tend to go into emission (locally) deep in the 


327 



atmosphere and therefore produce the consider- 
able backwarming observed. A rapid increase 
in r deep in the atmosphere increases the ten- 
dency of lines to go into emission and therefore 
increases the backwarming. 

The t term (Athay, 1972), which is propor- 
tional to Z, the net radiative bracket, tends to 
zero in the layers deep enough that the line is 
thermalized. Its principal effect appears near 
the surface. Although its value depends on the 
strength of the line and has somewhat different 
values for lines on the linear, flat, or damping 
part of the curve of growth, its value is certainly 
positive, so that it cools. For weak lines, c tends 
to exceed t and produces the cooling; for 
medium and strong lines, t exceeds c, often by 
a large amount, and is responsible for the sur- 
face cooling. Often, this cooling of outer layers 
is accomplished by a few very strong lines such 
as Na I D, Ca II H and K, Ca I 4227, Mg I b 
lines, and Mg II h and k lines; in other cases, 
the more numerous medium strong lines pro- 
vide most of the cooling except in the extreme 
outer atmosphere. 

What of the lines which warm the surface, 
such as TiO (Mould, 1975; Lengyel-Frey, 1977; 
Tsuji, 1978a; Krupp et al.; 1978)? Carbon 
(1979) explains these in a physical way based 
on analogy with processes in the continuum dis- 
cussed by Dumont and Heidmann (1973, 1976). 
In these processes, a major change in opacity 
occurs near the surface. Such processes can be 
discussed quantitatively within the above 
framework (Athay, 1972) if r o decreases rap- 
idly with optical depth. Some molecules (TiO, 
H 2 0, and HCN) exist in highest concentration 
near the outermost layers, where r o may 
change appreciably. In such cases, c becomes 
negative, and the surface layers are heated. This 
occurs for both the Schuster-Schwarzschild and 
Milne-Eddington atmospheres and is indepen- 
dent of the mechanism of line formation. The 
warming is enhanced by a flatter temperature 
gradient in the outer atmosphere, a line closer 
to LTE (e is larger), and a value of r Q which 
decreases with depth near the surface. An il- 


luminating example of an atmosphere with an 
opacity discontinuity has also been explored by 
Mihalas (1969). 

In actual calculation of stellar atmospheres, 
consideration of the detailed processes of line 
formation is normally avoided through the as- 
sumption of pure absorption (or even stronger 
assumption of LTE) so that the line source is 
equal to the Planck function for all lines, as 
shown in Figure 7-1. Occasionally, scattering 
is included for the strongest lines. Although this 
assumption may be adequate for deeper layers, 
its use might lead to fairly serious errors in the 
upper photosphere. Its only justification is a 
currently necessary simplicity. In fact, although 
some progress has been made, a thorough ex- 
amination of the mechanism of line formation 
and of the possible errors arising from the as- 
sumption of pure absorption is still lacking for 
even the strongest lines in the red-giant stars 
(Carbon, 1979; Gustafsson, 1981; de la Reza, 
this volume). Research on this problem can be 
expected to yield significant new knowledge. 

Treatments of Line Opacity 

Even with the assumption of pure absorp- 
tion in the line-source function, accounting for 
the rapid, almost random variation of the 
bound-bound absorption coefficient with wave- 
length across the entire spectrum, as implicit in 
the carbon-star spectrum shown in Figure 7-2 
(Querci and Querci, 1975b), is a formidable 
task, and some statistical procedure becomes 
necessary. The successive steps in treating such 
opacities form the backbone of a history of 
cool-star atmospheric modeling and the suc- 
cessful treatment of at least certain atomic and 
molecular opacities represents one of the signif- 
icant recent triumphs in stellar atmospheres. 

Intuitively, one knows that an increase of 
opacity in a given frequency band will inhibit 
the flow of radiation there, and flux constancy 
then requires an increase in radiative flux at 
neighboring frequencies. This leads naturally 
to the idea of a harmonic mean (HM) opacity 
(i.e., K~ l is the integral of the reciprocal of the 
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Figure 7-2: The observed spectrum of the N-type carbon star UU Aur in the region 6400 to 6480 cm 
(bottom graph) and the same spectrum corrected for weak telluric lines (middle graph). Many lines 
of C 2 and CN are identified (from Querci and Querci , 1975b). 


total absorption coefficient over a certain spec- 
tral interval divided by that interval). Unfortu- 
nately, different sources are not easily additive 
in an HM opacity, which makes its use cumber- 
some at best. Furthermore, the HM underesti- 
mates the opacity (Querci et al., 1972; Carbon, 
1974, 1979; see also Figure 7-4). It has been 
used in only one set of models (Auman, 1969), 
in which H 2 0 was the only opacity treated. 

A markedly different approach is repre- 
sented by the straight mean (SM) opacity, 
which is simply the average integrated absorp- 
tion coefficient of all lines in a spectral inter- 
val. An SM absorption coefficient (per mole- 
cule) for each species in each interval is fitted 


by a simple polynomial as a function of temper- 
ature, and the fitting coefficients are tabulated 
once for all (Alexander and Johnson, 1972; 
Johnson, 1974). Absorption coefficients for 
several absorbers are easily combined, and 
models are quickly calculated. The weakness of 
the SM is that the averaging process spreads the 
opacity over an entire spectral region, filling in 
the narrow spectral windows through which 
much of the radiation in reality escapes. The 
SM therefore overblankets the model, as has 
been elucidated by several studies (Querci et al., 
1972; Johnson, 1974; Carbon, 1974, 1979). The 
accuracy of the SM obviously improves as the 
spectral interval over which the smoothing is 
done is decreased. Johnson et al. (1972) show 
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an example for CN. For comparable spectral 
intervals, the SM can be thought of as a one- 
picket opacity distribution function (see also 
Figure 7-4). 

The opacity distribution function (ODF) 
was described by Strom and Kurucz (1966), and 
its modern use dates from about that time. In 
a given spectral interval (of, say, 50 A), one 
rearranges the absorption coefficient in wave- 
length to obtain a monotonic ordering, the 
result having the appearance of one giant spec- 
tral line (Figure 7-3). This new monotonic func- 
tion-termed the opacity distribution func- 
tion — is then approximated by a series of pick- 
ets as shown. (Care must be taken, of course, 
to flip the giant line from side to side in alter- 
nate intervals and depths to avoid systematic 
effects (Gustafsson et al., 1975).) Various math- 
ematical descriptions of the ODF are available 
(Querci et al., 1971; Carbon, 1973, 1979, 1984; 
Kurucz et al., 1974; Mihalas, 1978), and several 
applications have been made (Querci et al., 
1974; Querci and Querci, 1975a; Gustafsson et 
al., 1975; Kurucz, 1979; Eriksson et al., 1985). 


A representative ODF of CN and C 2 for 
conditions typical of carbon stars is shown in 
Figure 7-4 (Querci et al., 1971). Although the 
absorption of C 2 is much less than that of CN, 
it influences the ODF significantly. Also shown 
are the (constant) values for the SM and HM 
opacity; clearly the SM overestimates, and the 
HM underestimates, the true opacity. Similar 
comparison of SM, HM, and ODF opacities for 
CN in a carbon-rich model atmosphere was 
made by Carbon (1974). 

Calculating the ODF is tedious, but con- 
struction of models proceeds rapidly when the 
ODF is available. If, however, the ODF is sensi- 
tive in a given wavelength region to several 
opacity sources, the relative strengths of which 
might depend on a parameter (such as chemical 
composition) that is not known beforehand, 
one could be forced to a time-consuming itera- 
tive procedure (cf. Carbon, 1979, 1984). At 
least some of this liability could be overcome 
by a two-step procedure in which ODF’s are 
formed separately for each species and are then 



Figure 7-3. (a) Hypothetical spectrum over a spectral interval between X ; and \ r (b) Hypothetical 
line absorption coefficient over this interval, (c) The information in Figure 7-2b rearranged by strength 
of absorption to form an ODF or giant line (solid line) as described in the text. A possible six-picket 
representation of the giant line is shown by the dotted line, (d) The flux emitted by the picket represen- 
tation of the opacities in (c). 
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CUMULATIVE PROBABILITY ( P c ) 

Figure 7-4 . Comparison of various opacity 
treatments for CN and C 2 in the wave-number 
interval 6320 to 6420 cm" 1 . The graph shown 
(from Querci et al . , 1972) is for T = 2500 K, 
P = 1 dyne/cm 2 , and microturbulent veloci- 
ty VT = 5 km/s. 

reshuffled to account for the changing depth 
dependence of the opacities (Tsuji, 1976a). Re- 
cently, in what appears to be a significant ad- 
vance, the Swedish group has demonstrated the 
practicability of this suggestion by actual 
calculation and has further shown how the ac- 
curacy of the ODF can be increased by spline- 
fitting (Olander, 1981; Saxner and Gustafsson, 
1984). This method allows the ODF to accom- 
modate even the well-known problems of the 
depth-dependent combination of CN, C 2 , and 
CO opacities in cool carbon stars. 

Another drawback to the ODF arises from 
the rearranging of the absorption coefficient 
and its representation by a series of pickets; in 
this smoothing, all information regarding indi- 
vidual absorbers is lost. To ascertain the effect 
of a change in even one spectral line, as might 
arise from departures from LTE in even one 
atomic species, it is necessary not simply to re- 
compute an atmosphere with this changed 


value, but to recalculate the ODF. It is difficult 
even in principle to determine how this weak- 
ness of the ODF method can be overcome (Car- 
bon, 1984). For certain general problems, in 
which self-consistency between input parame- 
ters and results is necessary, the ODF may be 
too restrictive. For problems in which all pa- 
rameters to which the absorption is sensitive can 
be specified beforehand, however, the ODF 
method is convenient and accurate. In addition, 
the ODF appears to be capable of describing 
the opacity in a differentially moving atmo- 
sphere (cf. Mihalas et al., 1976). 

Suppose one begins to calculate model at- 
mospheres without averaging the opacity in any 
way (i.e., suppose one were not constrained by 
the amount of computing time available). Un- 
der these ideal circumstances, a straightforward 
procedure would be to calculate the integrated 
flux and its depth derivatives by a direct integra- 
tion over many monochromatic fluxes. Even 
though the variation of the absorption coeffi- 
cient with frequency is rapidly changing, its in- 
tegral can be obtained to any desired precision 
by taking a sufficient number of frequency 
points. This approach is called opacity sampling 
(OS), for there is no attempt to use all the 
opacity information; rather, one samples until 
one has sufficient. 

How many frequency points are needed? 
Very favorable answers are already available; 
even for complex opacities, only 500 to 1000 
frequency points are necessary if they are scat- 
tered over the entire spectrum and are randomly 
chosen with regard to opacity features (Peytre- 
mann, 1974; Johnson and Krupp, 1976). Mod- 
els computed with even as few as 100 frequency 
points match fairly well a more exact model in 
the photosphere, but depart from it in the outer 
layers, and several models with different sets 
of 100 frequencies lie within a certain envelope 
in a T-t diagram. As the number of frequency 
points is increased, the envelope narrows, and 
the envelope of models with 500 points has a 
spread which ranges from near zero deep in the 
photosphere to ±50 K near the surface 
(Johnson and Krupp, 1976). The final model 
is thus “converged” in the sense that, for more 
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than about 500 points, every set of points gives 
the same model. Related interesting numerical 
experiments with an ODF and an OS calcula- 
tion are reported by Gustafsson et al. (1975) 
and Gustafsson (private communication, 1984). 
Such characteristics put the OS method easily 
within the capabilities of current computers. 

The power of the OS method arises from its 
simplicity and its retention of all information 
on all lines which influence the transfer of ra- 
diation at the selected frequencies. Everything 
is under the control of the investigator. New 
values of such factors as chemical composition, 
turbulent velocity, or even departures from 
LTE are easily accommodated. This flexibility 
comes at a cost, however, because the construc- 
tion of each atmosphere demands much com- 
puting. Normally, the absorption of several tens 
of thousands of spectral lines (selected from a 
much larger sample) are accounted for in each 
calculation, making each iteration fairly time- 
consuming. Although some tricks (such as be- 
ginning with a small frequency set and using 
a larger set only for the final iterations) can save 
considerable computing, the OS method gener- 
ally requires more machine time per model than 
other methods. Furthermore, even after a 
model is converged, the computation of the 
emergent flux requires the calculation of fluxes 
at many frequencies within each element of 
spectral resolution, and the fluxes are then 
smoothed to the desired resolution. Finally, one 
must be cautious of OS models in the extreme 
outer photosphere, in which the few strongest 
spectral lines which control the energy balance 
and the thermal structure may not be properly 
sampled. Quantitative investigation of this 
point is needed. Use of the OS method for 
radiative transfer in an atmosphere with a ve- 
locity gradient is also very awkward (Carbon, 
1984) because of the shifting of line centers with 
depth. 

Although the ODF and OS may appear to 
be competitors, they are in some sense com- 
plementary (cf. Carbon, 1984). When all pa- 
rameters can be specified beforehand (e.g., in 
constructing a grid of stars of known composi- 
tion), the ODF method is preferable because of 


its speed. When absorption-affecting parame- 
ters are to be varied (e.g., in constructing a se- 
quence of stars differing in C/O or in account- 
ing for departures from LTE), the OS is su- 
perior. (Additional insight can be gained from 
reading Gustafsson et al., 1975; Johnson and 
Krupp, 1976; Carbon, 1979; Gustafsson, 1981; 
and especially Carbon, 1984; and Saxner and 
Gustafsson, 1984). We suspect that both the 
ODF and OS are capable of further develop- 
ment and look forward to these improvements. 

Less complicated than implied by the name 
(Voigt-analog/Elsasser-band model, or merci- 
fully VAEBM), an offspring of the ODF has 
been developed (Tsuji, 1976a) and used with 
success in constructing red-giant atmospheres. 
The method proceeds from a generalization 
(Golden, 1969) to Voigt profiles of the Elsasser 
band model, which treats the opacity of a com- 
plex molecule by a picket-fence model consist- 
ing of an infinite set of equally spaced spectral 
lines of equal strength. Although real spectral 
lines have different strengths, shapes, and 
widths, one can nevertheless define average val- 
ues of these quantities, adding pickets as 
needed. Even overlapping lines can be roughly 
accounted for. Restrictive as these assumptions 
appear, the VAEBM representation is, over a 
spectral interval sufficiently small that the 
Planck function does not vary appreciably, for- 
mally equivalent to the ODF described earlier 
(Tsuji, 1976a), and it has been exploited with 
good results. However, one can easily imagine 
the method failing in the case of overlapping 
opacities, especially with different depth de- 
pendences (Saxner and Gustafsson, 1984). In 
short, the VAEBM is a simple ODF, with the 
strengths and weaknesses of the ODF in addi- 
tion to the more severe approximations by 
which the VAEBM gains its simplicity. On the 
other hand, it has the advantages of speed and 
the power to treat, at least crudely, even such 
complicated molecules as H 2 0. Certain 
molecular VAEBM opacities have been com- 
pared to Rosseland (harmonic mean) opacities 
and JOLA opacities (just overlapping line ap- 
proximation) in an interesting paper (Tsuji, 
1971). 
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Like a phoenix, SM opacities never really 
die, and they have recently been in the news 
again because of their use in spherical models. 
A comparison of Rosseland opacities from SM 
and OS treatments (Scholz and Tsuji, 1984) 
finds good agreement for solar composition, 
but for a carbon-rich composition, the OS 
opacities are much smaller. Unfortunately, the 
agreement for solar composition is largely for- 
tuitous because the OS treatment included 
H 2 0 as SM opacity (Johnson et al., 1980), and 
at the temperature used ( T = 3000 K), the OS 
opacities are in fact dominated by H 2 0. Ac- 
cord between “OS” and “SM” is therefore 
simply accord between the SM of Johnson and 
the SM of Tsuji. Clearly, the Rosseland opacity 
must always be greater for SM than for OS or 
ODF, because the SM approximation fills in all 
opacity windows and acts, in a given spectral 
interval, exactly as an additional continuous 
opacity source. 

In the following thought experiment, we 
imagine a continuous opacity and a superposed 
line opacity. Let the lines (rectangular for sim- 
plicity) cover half the spectrum and have an 
(equal) absorption coefficient 100 times that of 
the continuum. If calculated by the ODF or OS 
method, the Rosseland opacity is approxi- 
mately twice that of the continuum alone since 
that half of the interval covered by lines is ex- 
tremely opaque. If calculated by the SM, the 
average line opacity would be 50 times the con- 
tinuous opacity, and the Rosseland opacity 
would also be increased by a factor of 50 . 
Clearly, Rosseland opacities and SM opacities 
are not to be mixed . Yet consider the absorp- 
tion of a molecule with enormous numbers of 
weak lines, such as H 2 0 or HCN. If the 
100000 strongest lines of HCN make little dif- 
ference in an atmosphere (Tsuji, 1984), but 6 
million weak lines make a large difference 
(J0rgensen, 1985), we understand the reason. 
Such an enormous number of lines fills in all 
the spectral windows, and like the SM, acts at 
that depth as an additional continuous opaci- 
ty. If a molecule contributes a haze of lines 
(even weak ones), it may greatly increase the 
opacity, and such polyatomic opacities 


(Auman, 1969; Tsuji, 1984; Eriksson et al., 
1984) must be included (laboriously) in at least 
the coolest red-giant atmospheres. A similarly 
large Rosseland opacity can be obtained by fill- 
ing in the opacity windows with isotopic lines 
and very large microturbulant broadening 
(Sharp, private communication, 1984). 

THEORETICAL 
PHOTOSPHERIC MODELS 

As might be expected, the passage of time 
has seen a steady improvement in theoretical 
photospheric models, particularly in regard to 
the treatment of opacity. It is convenient, then, 
to separate models into two groups according 
to the opacity approximation employed: (1) 
first-approximation opacities (SM, HM, or sim- 
ilar treatment), and (2) second-approximation 
opacities (ODF, OS, or VAEBM treatments). 
We include in this chapter only giant and 
supergiant models with effective temperatures 
below 4000 K. 

First-approximation models have already 
been summarized and extensively discussed 
(Vardya, 1970; Johnson, 1972, 1974; Kipper, 
1973; Carbon, 1979; Gustafsson, 1981). 
Although they are now generally superseded by 
superior models with the same range of 
parameters, these earlier models were the basis 
of several noteworthy discoveries. Perhaps the 
most significant was that, contrary to common 
belief, convection (in the usual mixing-length 
formulation) was of minor importance in the 
photosphere itself— both for M stars (Auman, 
1969) and for carbon stars (Johnson, 1972). 
Somewhat surprisingly, even though models 
with SM and HM opacities have quite different 
structures, computed molecular column den- 
sities agree very well (Goon and Auman, 1970; 
Johnson et al., 1975), indicating that these lat- 
ter are rather insensitive to the details of the 
atmospheric structure. The large scale heights 
observed in stars such as f Aur were duplicated 
by theoretical models with appropriately low 
values of surface gravity (Johnson, 1972). The 
observed flux from a cool carbon star was fairly 
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well represented by the predicted flux from an 
appropriate model (Johnson, 1972), except that 
the model flux was too high in the blue and the 
CN bands were too deep. Older atmospheric 
models were used for several projects: (1) to ex- 
amine the remarkable cooling by CO in the 
outer layers of stars (Johnson, 1973); (2) to in- 
terpret the infrared spectra of several carbon 
Mira stars in support of the hypothesis of C 
enrichment rather than CNO processing 
(Thompson, 1977); (3) to study the atmosphere 
of Betelgeuse (Fdy and Johnson, 1973); (4) to 
deduce, from IR observations, approximate ef- 
fective temperatures of K and M giants and 
supergiants (Scargle and Strecker, 1979); and 
(5) to infer the ratio 12 C/ 13 C in five carbon 
stars (Olson and Richer, 1979). Finally, to ex- 
plore the gross molecular features of S stars, 
Piccirillo (1980) computed 85 SM model at- 
mospheres for 2500 = T eff ^ 3500 K and log 
g = 0 for compositions enriched in C to give 
ratios of C/O in the range 0.6 (solar) to 1.00. 
A point of interest here is his finding that the 
thermal structures of some of these SM mod- 
els were similar to OS models. Several conclu- 
sions regarding C/O and s-process enhance- 
ments in S stars emerged from a comparison 
of his molecular column densities and observed 
band strengths — perhaps the most significant 
being that, in certain S stars, C/O may be 
nearer 0.9 than 1.0. 

In considering models with second-genera- 
tion opacities , we continue to distinguish 
models which are oxygen-rich (M stars) and car- 
bon-rich (C-stars). S stars, which are intermedi- 
ate, will be discussed along with C stars. 

Models for M Stars 

An overview of currently available models 
for solar composition is given in Figure 7-5 
(from Johnson, 1985), where models are shown 
in the parameter space of effective temperature 
and surface gravity (log g). All published 
models with the assumptions of Figure 7-1 and 
plane-parallel geometry (PPG) and based on 
ODF, VAEBM, or OS opacities are shown; 
they are identified as BEGN (Gustafsson et al., 
1975; Bell et al., 1976a). T (Tsuji, 1976a; 
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Figure 7-5. Overview of available models for 
normal red giants of solar composition dis- 
played in the space defined by the surface grav- 
ity ( log g) and the effective temperature. Refer- 
ences are given in the text (from Johnson, 
1985). 

1978a), or JBK (Johnson et al., 1980). In each 
group, certain additional models with varia- 
tions in parameters other than those shown as 
coordinates were calculated. Models of Tsuji 
(1976a, 1978a) are purely radiative, whereas 
models of Gustafsson et al. (1975), Bell et al. 
(1976a), and Johnson et al. (1980) also include 
convection. Although it may appear that this 
parameter space is fairly well covered by 
models, the coverage is sparse in places, and 
several principles of the classic atmosphere may 
be stretched rather badly toward lower 
temperatures. Perhaps the worst assumption is 
the crude treatment of H z O opacity, but the 
assumption of PPG and even homogeneity may 
not be met in the coolest, most extended stars. 

Designed to represent Population I and II 
giants, the grid of Gustafsson et al. (1975) and 
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Bell et al. (1976a) spans the ranges 3750 K S 
T ff g 6000, 0.75 < log g < 3.0, and -3.0 5 
[A/H] s 0.0. Altogether, 103 models based on 
ODF opacities for CO, CN, OH, NH, MgH, 
and atomic lines were computed. No TiO lines 
were included, which unfortunately compro- 
mises the accuracy of models with r ff S 4000 
K — precisely those of interest here. Convection 
was treated by the mixing-length theory with 
a mixing length 1.5 times the pressure scale 
height. Although turbulent pressure was ne- 
glected, a study of the possible effects of its 
neglect was made. For a model from near the 
middle of their grid, Gustafsson et al. (1975) 
find that, for u(turb) = 5 km/s, P(turb) 
becomes an important component of P(total), 
but the thermal structure [P(r)] is affected on- 
ly slightly except in the convection zone. The 
surface cooling of CO lines was confirmed 
quantitatively, and CN lines were shown always 
to backwarm the entire atmosphere. The 
authors conclude (Gustafsson et al., 1975) that 
self-consistent LTE models with as many as 
10 5 wavelength points are possible, but a much 
smaller number of points is satisfactory for 
either the ODF or OS method. The temperature 
uncertainty due to all accumulated errors is 
estimated at 100 K in the outer layers, 20 to 30 
K through most of the photosphere, and 100 
K in the convective zone. Uncertainties due to 
the assumptions of LTE and homogeneity, 
which were not investigated, could be con- 
siderably larger. Several applications of these 
models to astrophysical questions are described 
in the section Spectral Distribution of Energy 
Fluxes (see also Gustafsson, 1981). 

To find a middle ground between the crude 
approximations inherent in SM opacities and 
the tedious computations required for the usual 
ODF opacities, Tsuji (1976a) used the VAEBM 
method described earlier to treat the opacities 
of OH, CN, CH, TiO, MgH, SiH, CaH, CO, 
and H 2 0 (some approximated more carefully 
than others) and construct 11 models for 
supergiant stars. They range in effective 
temperature from 3400 to 4000 K and log g = 
0.0 or -0.5. Turbulent pressure is included in 
the hydrostatic equilibrium equation with 


velocities of 3, 6, and 12 km/s, but convection 
is neglected. In a companion paper, Tsuji 
(1976b) uses these models for a specific discus- 
sion of Betelgeuse, for which he deduces a value 
of 7\ = 3900 ± 150 K. 

eff 

A grid of 22 models for K and M giants 
based on these VAEBM opacities was later 
computed and used to define a new relation be- 
tween spectral type and effective temperature 
(Tsuji, 1978a). Based on LTE, hydrostatic equi- 
librium (including radiative and turbulent pres- 
sure with a turbulent velocity of 3 km/s), and 
radiative equilibrium (convection was ne- 
glected), these models were computed for the 
ranges 2600 ^ T ff S 4200 K and -2.0 § log 
g < 2.5. Predicted fluxes were then compared 
to photometry of a number of cool stars to fix 
the relation between effective temperature 
(from the models) and spectral class (from the 
stars). Tsuji emphasizes the desirability of com- 
paring infrared fluxes for finding best fits be- 
cause of the relative insensitivity of the infrared 
flux to opacity variations. Yet we wonder 
whether a fit also to the visual region (with its 
much steeper slope) as well might fix the 
temperature of a star with greater precision. 
Many of these models predict substantially too 
much flux in the visual region (a common 
problem). 

Based on the ATLAS computer program 
(Kurucz, 1970), Johnson et al. (1980) calculate 
40 atmospheres with solar composition and the 
parameters: T cff = 2500, 2750, 3000, 3200, 
3400, 3600, 3800, and 4000 K and log g = 0.0, 
0.5, 1.0, 1.5, and 2.0. Convection is included, 
but turbulent pressure is neglected. Atomic lines 
and molecular lines from CN, CO, C 2 , TiO, 
OH, NH, CH, and MgH are treated by the OS 
method; H 2 0 is added as an SM opacity. In a 
sense, these results complement those of Gus- 
tafsson et al. (1975) to provide a complete set 
of models for cool giants and supergiants for 
all temperatures. Although this work is the 
most complete now available for very cool 
stars, the weakness of several assumptions 
should caution the reader against uncritical ac- 
ceptance of every value. Treating H 2 0 as an 
SM opacity overblankets the coolest models 


335 



(r eff < 3200 K) (Auman, 1969; Tsuji, 1971), 
and this is probably the most serious defect. 
Neglect of grain formation in the coolest 
models may also be incorrect. The temperature 
and density in the surface layers of the super- 
giant (log g = 0.0) models might be even lower 
if sphericity were included. 

How well do the models agree where they 
overlap? A comparison of T(j ) has been made, 
for the particular case of (4200/2. 25/solar), be- 
tween a Bell et al. (1976a) (ODF) model and 
a Johnson et al. (1980) model (Johnson and 
Krupp, 1976). All opacity parameters were 
identical, but the opacity treatment, as well as 
the computer codes, was different. The ex- 
cellent accord in T(f) — within ± 30 K except 
in the extreme outer layers — indicates that 
both the ODF and OS methods are satisfactory. 
In the extreme outermost layers, the tempera- 
ture is fixed by those few strong lines which still 
have some opacity (i.e., for which t(v) * 1), 
and these may not be accurately treated by the 
OS method (Carbon, 1979; Johnson et al., 
1980). Departures from LTE and possible chro- 
mospheric heating are also more likely in these 
superficial regions. Since TiO opacity, which 
is important in warming the outer layers, is ne- 
glected by Gustafsson et al. (1975) and Bell et 
al. (1976a), their coolest models ( T ff < 4000 
K) are not recommended (McGregor, 1980). 

A parallel comparison of a model (4000/ 
1.5/S) with VAEBM opacities from Tsuji’s grid 
with a similar model of Bell et al. (1976a) shows 
that the two agree within ± 50 over most of the 
atmosphere, but the former is hotter in the 
outer layers and cooler in the deeper layers (by 
up to 130 K). Emergent fluxes from the two 
also agree well, whereas similar fluxes from an 
SM model display molecular bands which are 
too deep (Tsuji, 1976a). 

The atmospheres of Johnson et al. (1980) 
form a sequence parallel to that of Tsuji 
(1978a), and together these constitute the only 
models for M giant stars with realistic opacities. 
A comparison of the T-t relation between 
models from each set has been made by John- 
son et al. for the parameters 3600/0.0/S, and 
the agreement is excellent. This is the more im- 


pressive because the computer codes and molec- 
ular data were independent and the methods of 
treating the line opacities were different 
(VAEBM by Tsuji and OS by Johnson et al.). 

Models for Carbon and S Stars 

An overview of the available models for car- 
bon and S stars (i.e., with C/O > 0.6) is pro- 
vided by Figure 7-6 (Johnson, 1985), which 
shows the regions in the parameter space of ef- 
fective temperature and the ratio of carbon to 
oxygen covered by various sets of models. Fig- 
ure 7-6 is, of course, a two-dimensional projec- 
tion of the multidimensional space whose axes 
are the parameters specifying these models. Be- 
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Figure 7-6. Overview of available models for 
peculiar red giants (M, S, R, and N-type stars) 
displayed in the space defined by the carbon- 
to-oxygen ratio and the effective temperature. 
References are given in the text (from Johnson, 
1985). 
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sides surface gravity, they include the abun- 
dance of all significant elements, especially C, 
N, and O. Clearly, Figure 7-6 is a considerable 
simplification, and sets of models shown as 
lines (such as QQ) contain models varying in 
some parameter other than the two used as co- 
ordinates— in this case, N/H and log g. 

Based on an ODF which included lines of 
the molecules CN, C 2 , and CO, but no atomic 
lines, the French group (Querci et al., 1974) 
constructed 12 carbon-star model atmospheres 
with effective temperatures of 4500, 4200, 3800, 
and 3400 K and surface gravities of g = 0.1, 
1.0, 10 cm/s 2 . These were later supplemented 
with 9 cooler models having effective tempera- 
tures of 3000, 2600, and 2200 K (Querci and 
Querci, 1974). Convection was neglected, but 
a turbulent pressure P(turb) = /3q v 2 was used 
in the hydrostatic equation. Values of f3 = 1 
and (3 = 0 and u = 5 km/s were chosen to agree 
with the microturbulent velocity deduced earlier 
from the spectrum of UU Aur. Solar abun- 
dances were used for all elements except CNO, 
which were fixed at H/C/N/O = 1/4.1 x 
10 _5 /1 .48 x 1 0' 3 / 1.25 x 1 0' 5 (C/O = 3.2 and 
N/C = 37), to represent a composition pro- 
duced by CNO processing and mixing. The ef- 
fect of molecule line blanketing is clearly seen 
in their models (Figure 7-7); for successively 
lower temperatures, the molecular line blanket- 
ing increasingly cools the surface layers and 
warms the photosphere compared to the un- 
blanketed models. These authors found: (1) 
considerable surface cooling and photospheric 
heating due to CO, C 2 , and CN; (2) good 
agreement in thermal structure (perhaps sur- 
prising) with certain SM models; (3) fair agree- 
ment— slightly better than with a blackbody— 
of predicted fluxes, except in the visual region, 
with wideband photometry of TX Psc for mod- 
els with J ff = 3000 and 3400 K. An addi- 
tional set of 14 models with compositions rep- 
resenting both CNO cycling and triple-alpha 
processing were calculated, and certain molec- 
ular features characterizing each composition 
were noted (Querci and Querci, 1975a). Syn- 
thetic spectra from these models are discussed 
later. 



log T s td 

Figure 7-7. Effect of molecular line opacity. At 
an effective temperature of 3400 K and surface 
gravity of 1.0, a photospheric model for a car- 
bon star with C 2 , CN, and CO opacity is com- 
pared to an identical model without these 
(Querci and Querci, 1974). In each case, the 
blanketed model is cooler at the surface and 
warmer in the photosphere. 


Another set of 27 models for M, S, and C 
stars has been produced with OS opacities for 
atomic lines and the molecules CH, NH, OH, 
MgH, CO, CN, and C 2 (Johnson, 1982). 
Three series of models were computed, all with 
log g = 0.0, at effective temperatures of 2500, 
3000, and 3500 K. Beginning with a model of 
solar composition, carbon is added incremen- 
tally to give C/O ratios of 0.60 (approximately 
solar), 0.90, 0.95, 0.98, 1.00, 1.02, 1.05, 1.20, 
and 2.00. Convection is included in the usual 
mixing length formulation with 2 / H - 1.0; 
turbulent pressure is neglected. The effect of 
C/O on the structure is illustrated in Figure 7-8, 
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which displays models with C/O = 0.60, 1.00, 
and 2.00 — possibly representing a cool M su- 
pergiant, an S star, and an N-type carbon star, 
respectively. For comparison, a model with so- 
lar composition but no line blanketing (labeled 
0.60*) has the expected form: a steep gradient 
in the photosphere and a flat outer portion with 
a boundary temperature T q = 2300 K. An 
identical model that includes atomic and molec- 
ular opacities is warmer throughout the photo- 
sphere due to the backwarming and cooler (due 
mainly to CO and H z O) in the outer layers 
(outside the graph shown) with T q = 2100 K. 
As carbon is added and free oxygen is reduced 
by the formation of CO, the partial pressures 
and opacities of H z O and TiO decrease. Recall 
that H 2 0 cools near the surface and back- 
warms the photosphere, while TiO warms the 
surface layers (Krupp et al., 1978; Gustafsson 
and Olander, 1979; Gustafsson, 1981). This de- 
crease in backwarming as H z O and TiO 
decrease in abundance causes the temperature 
drop noted. In the outer atmosphere, the in- 
crease in CO cooling as carbon is added more 
than compensates for the loss of H z O cooling, 
and the outer layers continue to cool (Carbon, 
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Figure 7-8. Effect of the enhancement of C on 
the thermal structure (T(t / )) of a model with 
(3000/0.0/solar except C) (from Johnson, 
1982). Optical depth is measured in the contin- 
uum at 1 \n.m. The parameter is C/ O. A dotted 
line is a model with solar composition (C/O = 
0.60) but without any line blanketing. 


1974; Heasley et al., 1978; Carbon, 1979). The 
remarkable resemblance between the model 
with C/O = 1.0 and the unblanketed model il- 
lustrates well that an S star has little molecular 
line opacity. 

As additional carbon is added beyond C/O 
= 1.0, the temperature increases again 

throughout the photosphere due to the back- 
warming of CN and C 2 (Querci et al., 1972, 
1974; Querci and Querci, 1975a; Sneden et al., 
1976). In the outer layers, the temperature 
seems to decrease. In actuality, the boundary 
temperature (1400 K) does not fall, but it is 
reached at a larger value of optical depth as 
C/O increases; that is, this cool outer layer 
deepens as carbon is enhanced (Johnson, 1982). 
Figure 7-8 illustrates that, because an enhance- 
ment of carbon changes the atmospheric struc- 
ture, spectral features — including those of 
elements other than carbon or its compounds — 
may also change. This point has previously 
been made from a study of molecular column 
densities (Johnson et al., 1975; Piccirillo, 1980; 
Johnson and Sauval, 1982), approximate band 
strengths (Scalo, 1973), and synthetic spectra 
(Querci and Querci, 1976). 

Of great interest is the value of the boundary 
temperature and the thermal structure of the 
outer layers of red-giant atmospheres. In his re- 
view of model atmospheres for intermediate 
and late-type stars, Carbon (1979) pointed out 
that, in the discovery paper (Johnson, 1973), 
the cooling due to CO (based on SM models) 
was twice that found in later studies with more 
accurate ODF opacities (Gustafsson et al., 
1975). In this present case, the cooling of the 
M supergiant is about 200 K, whereas that of 
the S star and carbon star are about 900 K. In 
addition, Johnson et al. (1980) list three consid- 
erations which might further lower the bound- 
ary temperature: (1) models may not be calcu- 
lated at sufficiently small optical depths to 
reach the true value of T o ; (2) the OS method 
may fail to sample the strongest lines in the out- 
ermost layers and therefore underestimate the 
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cooling; and (3) neglected sphericity effects may 
lower T q . A cooler surface layer will increase 
the chance for grain formation in the photo- 
sphere (Lucy, 1976; Menietti and Fix, 1978; 
Deguchi, 1980; Schmid-Burgk and Scholz, 
1981; Draine, 1981; Alexander et al., 1983). On 
the other hand, chromospheric heating may 
raise the boundary temperature. At present, the 
boundary temperatures for these stars are un- 
known, but the value of T . may be obtain- 
able from high-resolution International Ultra- 
violet Explorer (IUE) observations, and such 
investigations are strongly encouraged. 

Although the R stars, at least those of spec- 
tral class R0 to R3, correspond more closely to 
K than M stars, they are included in this vol- 
ume because of their carbon-rich spectra and 
advanced evolutionary state. Models have been 
computed by two groups. The Swedish group 
published 20 models (Olander, 1981) with effec- 
tive temperatures in the range 3800 to 4800 K. 
An interesting aspect of this project is the at- 
tempt to circumvent some of the drawbacks of 
the ODF (Carbon, 1979) by reshuffling the 
ODF’s of different species. The author con- 
cludes that: (1) because of the relatively high 
temperatures of R stars, changes in composi- 
tions do not drastically affect the temperature 
structure; (2) raising the N abundance causes 
a slight backwarming (< 100 K) and surface 
cooling due to the increased CN opacity; (3) the 
effect of changes in the O abundance is very 
slight; and (4) raising the microturbulent veloc- 
ity from 2 to 5 km/s causes a slight warming 
throughout the atmosphere, as was found for 
G and K giants by Gustafsson et al. (1975) and 
Bell et al. (1976a) and for cooler carbon stars 
by Johnson (1982). At approximately the same 
time, the Indiana group produced 1 1 models for 
R stars (Johnson and Yorka, 1985). These span 
a slightly higher temperature range (4200 to 
5600 K), and surface gravity is fixed at either 
log g = 2.0 or 3.0. The log abundances of 
H/C/N/O were taken as 12.00/9.09/8.31/8.55, 
to give C/O = 1.74, although models were 
calculated for the additional values 0.60, 1.00, 
and 2.50 at a temperature of 4600 K. Yorka 


(1981) found the two sets of models to be in 
excellent agreement, and Dominy (1984) used 
them to help infer the effective temperatures of 
a set of warm R stars (R0 to R3). 

None of the carbon-star models published 
to date contains polyatomic opacities (except 
for a few SM models by Johnson, 1974), but 
the importance of HCN in cool carbon stars has 
recently been reemphasized (Eriksson et al., 
1984; Tsuji, 1984; J0rgensen, 1984). A large set 
of models for 2500 S P ff 5 3500 K and 1.00 
S C/O £ 50 which include ODF opacities of 
HCN and C 2 H 2 , as well as the usual diatomic 
opacities, has been calculated. We refer to these 
as the EGJN models (Ericksson et al., 1985) 
and eagerly await their publication. 

Finally, we mention a set of 10 carbon-rich 
models with the parameters 3000/0.0/S except 
C; C/O = 1.05, in which the fractional abun- 
dance of H was incrementally reduced and re- 
placed by He (Johnson et al., 1985). Interest- 
ingly, even this drastic change in composition 
alters the thermal structure ( T-t ) only very 
slightly, whereas the pressure structure (P-7) is 
greatly changed (equivalent to a change in grav- 
ity). The energy flux curve is also changed in 
the sense that, as the hydrogen abundance is 
decreased, the H~ flux peak at 1.67 fim is 
steadily reduced relative to the blackbody peak 
near 1.0 fim. 

Again, how well do the model makers agree 
among themselves? A comparison of models 
for carbon stars is more difficult than for M 
stars because the extra dimensions of chemical 
composition mean less overlap between differ- 
ent sets of models. Nevertheless, in Figure 7-9 
we compare, for models with the parameters 
3000/0.0/solar except C; C/O = 2.00, the ther- 
mal structure ( T(t )) and the pressure structure 
( T(P )) of the HJ (Johnson, 1982) models and 
EGJN models (Gustafsson, private communi- 
cation, 1984). For models with diatomic mo- 
lecular opacities only, both the thermal struc- 
ture and the pressure structure are in good 
agreement, the differences being attributable to 
differences in molecular data and computer 
codes. When polyatomic opacities are includ- 
ed, the thermal structure is almost unchanged, 
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Figure 7-9. Comparison of carbon-rich models computed by the Swedish group (EGJN) and the Indi- 
ana group (HJ) and the effect of polyatomic opacities . (a) The T-t structure of models by the two 
groups with (EGJN poly) and without (HJ, EGJN) HCN opacity, (b) The T-P structure of models 
with and without HCN opacity. ( Gustafsson , private communiction , 1984; other references are given 
in the text.) 


but the pressure is greatly reduced because of 
the large increase in opacity. 

From this brief overview of models with 
carbon-rich composition, it should be apparent 
that the coverage of parameter space is sparce; 
we are barely beginning. Before using any 
model atmosphere for red-giant stars to calcu- 
late quantities, such as profiles of very strong 
lines, which depend sensitively on the tempera- 
tures in the extreme outer layers (r < 0.0001 , 
say), one should contact the maker of the model 
to discuss problems. (The guarantee normally 
covers only the photosphere!) 


SPHERICAL MODELS 

While calculations (Auman, 1969; Johnson, 
1972; Tsuji, 1976a) confirm the adequacy of 
PPG for giants (log g = 2.0), it may fail for 
low-mass supergiants and Miras. Three changes 
must be made in the usual equations to accom- 
modate sphericity. 

1. The equation of radiative transfer must 

include a curvature term: 


9 / (1 -n 1 2 ) 9 / 


where /(/*,/*,!') is the specific intensity, g, 
= cos 0 as usual, e(r,^) is the volume 
emissivity, and x(r 9 v) is the linear extinc- 
tion coefficient. Azimuthal symmetry 
and time-independence are assumed. If 
(dl/d{i) = 0, one recovers the usual equa- 
tion of PPG. 

2. The surface gravity is no longer constant, 
but g - GM/r 2 , and the variation of g 
(but not M) with r must be accounted 
for. 

3. A third parameter, in addition to T ff 
and g, is needed to specify fully the atmo- 
sphere. This is often the stellar radius, 
measured to some reference level, such 
as t = 1.00. 

The complete equations for spherical 
geometry and methods for their solution have 
been available for several years (cf. Mihalas, 
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1978). Rather different mechanisms produce at- 
mospheric extension in the early-type and the 
late-type stars. In hot stars, electron scattering 
produces an outward radiative force which 
leads to an extended atmosphere. Contrarily, 
red-giant atmospheres are extended only be- 
cause of the very low surface gravity (unless, 
of course, there are outward flows; see Gold- 
berg, this volume). 

A commonly used measure of the impor- 
tance of sphericity is the extension , or the 
relative thickness of the atmosphere compared 
to the stellar radius. More precisely, the exten- 
sion is d = (/*( top) - r(phot))/r(phot), where 
r(top) is the radius measured to the point at 
which t - 0.0001, say, and r (phot) is the stellar 
radius measured to r = 1.00. Authors differ 
slightly in employing a Rosseland optical depth 
or a monochromatic optical depth and in the 
exact value of the optical depth (10~ 5 to 10" 3 ) 
chosen for the top of the atmosphere. 

To delineate the regions of the HR diagram 
in which sphericity might be important, Hundt 
et al. (1975) and Schmid-Burgk and Scholz 
(1975) computed a series of spherical models 
for late-type stars having gray opacities. Of in- 
terest here are two principal findings: (1) many 
red giants and supergiants of moderate mass ( M 
~ 1 A^7 o ) develop extended atmospheres at 
some epoch during their lives, usually near the 
peak of their RG luminosity; and (2) tempera- 
tures in the outer layers are considerably lower 
than in a comparable plane-parallel model. 

More realistic models were obtained by Wat- 
anabe and Kodaira (1978, 1979), who included 
the line opacities of H 2 0, CO, CN, OH — ap- 
proximated with a mixture of SM, JO A (just 
overlapping approximation), and VAEBM. In 
these spherical models, the temperature was 
lower in the outer layers (by 150 to 200 K in 
the model with T ff = 3200 K, for example) 
and higher in the deeper layers than in a corre- 
sponding PPG (here called compact) model. 
The density was also lower in the outer layers 
than in a corresponding PPG model. The num- 
ber density of most molecular species, such as 
CO and CN, followed the overall decrease, but 
there were certain exceptions, notably those for 


which the increased efficiency of molecular for- 
mation due to the lower temperature overcomes 
the reduction in the overall density (e.g., H 2 0 
in the model with 7 ff = 3200 K). For at- 
mospheres with extensions up to d = 1.5, these 
workers found: (1) an increase of atmospheric 
extent, (2) a reddening of (R-I), and (3) a 
decrease in the absolute visual magnitude. In 
addition, stellar diameters measured at 
wavelengths within those molecular bands 
formed in extreme outer layers and hence sen- 
sitive to atmospheric extension, such as H 2 0 
and TiO, increased as much as a factor of 1.5 
compared to the adjacent continuum. 

Similar calculations which reinforced these 
tentative conclusions were carried out by 
Schmid-Burgk et al. (1981) and Kipper (1982) 
for models with solar composition. Most of the 
conclusions were drawn from models computed 
with continuous opacities and H 2 0 only, 
which is probably the dominant absorber in this 
temperature range. These authors document the 
great sensitivity of the thermal structure to the 
treatment of molecular opacities, the incorrect 
treatment of which can obscure other effects. 
Figure 7-10 (Schmid-Burgk et al., 1981) il- 
lustrates one of the principal effects of 
sphericity — a lowering of the temperature in the 
outer layers, where divergence of the flux is im- 
portant. The temperature drop is greater in the 
lower gravity A models. (A slight backwarm- 
ing — too small to be seen here— accompanies 
the surface cooling.) In the outermost layers, 
the temperature decrease is 300 to 500 K for a 
model with extension d = 0.25 compared to a 
compact model. The other principal effect of 
atmospheric extension— the decrease in density 
in the outer layers— is illustrated in Figure 7-1 1 
(Schmid-Burgk et al., 1981). The key effect is 
the drop in pressure (and density as well) from 
A0 to A2 and B0 to B3; the lower gas density 
arises from the decreasing gravity, which falls 
as r" 2 in these spherical models. Figure 7-11 
also shows the radius (r/R ) plotted against an 
optical depth scale (7(1.2 ^m)); the trans- 
formation from radius to optical depth is sensi- 
tive to the details of the opacities and to the 
relation of H and H 2 , which dominate the 
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Figure 7-10. The effect of atmospheric exten- 
sion on the thermal structure of a star. For each 
series of models, the extension of the atmo- 
sphere is increased from AO (PPG) and BO 
(PPG) to A3 (d = 1.39) and B3 (d = 0.20). 
Parameters of the models are: A: 3 000/ -1.0; 
B: 3000/0.0 (from Schmid-Burgk et at, 1981). 

equation of state at these low temperatures 
(Schmid-Burgk et al., 1981). As is clear from 
a consideration of the two figures, the number 
density or partial pressure of a particular spe- 
cies depends rather critically on the details of 
the temperature/density relation in the outer at- 
mosphere, which is not well known at this stage 
of research. Schmid-Burgk et al. also estimate 
the maximum tolerable value of the extension 
to be d = 0.05 or perhaps as large as d = 0.10 
in certain cases. The values are exceeded, at 
some epoch, in almost all M giants with lumi- 
nosities g 10 3 solar luminosities and masses 
near 1 M© (see, for example, Scalo et al. 
(1978) and Scalo and Miller (1979) for an 
estimate of these masses). 



Figure 7-11. Effect of atmospheric extension on 
the gas pressure for the models of Figure 7-10 
(from Schmid-Burgk et al., 1981). The radius 
(r/R) is also shown for the most extended 
models. 

Pursuing further the consequences of varia- 
tions in the chemical composition for M super- 
giants, still within the framework of SM molec- 
ular opacities, Wehrse (1981) varied separately 
the abundances of He, the CNO group, and the 
metals. Lowering the metal abundance in- 
creases atmospheric extension and surface cool- 
ing, whereas introducing a hydrogen deficiency 
reduces sphericity effects. However, because of 
its interplay of various factors, separating 
sphericity from compositional effects appears 
to be extremely difficult. 

To find an observational third parameter, 
Scholz and Wehrse (1982) computed very cool 
(T ff = 2750, 3000 K) models based on SM 
opacities of several molecules in addition to 
H 2 0 and then calculated band strengths from 
these models. They find that for their models : 
(1) the 1.04-4.00 and the 2.10-4.00 /tm colors 
can be calibrated to yield effective tempera- 
tures; (2) certain metal/ion atomic lines can 
yield the surface gravity; and (3) TiO lines in 
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the red are especially sensitive to atmospheric 
extension. 

Again utilizing SM opacities, Scholz and 
Tsuji (1984) study carbon-rich spherical atmo- 
sphere. Because of the greater temperature sen- 
sitivity of the carbon-rich opacities, these atmo- 
spheres are less extended than oxygen-rich at- 
mospheres. For the same reason, gas pressure 
and density are higher in C-rich spherical atmo- 
spheres than comparable PPG (compact) atmo- 
spheres, but oxygen-rich stars behave 
oppositely. 

Spherical models will clearly come to play 
a major role in the study of red-giant stars — 
especially for supergiants and Miras — as soon 
as realistic opacities can be used. 

CONVECTION AND 

INHOMOGENEITIES 

Despite much research on the subject, con- 
vection is one of the more uncertain aspects of 
atmospheric modeling. The problem is simply 
that we do not understand the phenomenon; 
not only is convection an extremely complex hy- 
drodynamic phenomenon, but we are without 
the benefit of sure terrestrial or laboratory 
analogs. 

Envelopes of red-giant stars must be the 
scene of deep strong convection due to the 
opacity associated with hydrogen and helium 
ionization/recombination zones. Somewhere in 
the deep photosphere — near the depth where 
t(R) = 10, but depending on effective tempera- 
ture, surface gravity, and composition — the en- 
velope changes rapidly from a convective to a 
radiative regime. Photospheres lie at and imme- 
diately above this transition, and it has been 
shown for both M stars (Auman, 1969) and car- 
bon stars (Johnson, 1972) that convection, if 
treated in the usual local mixing-length theory 
without overshoot, is unimportant in the photo- 
sphere (see also Gustafsson et al., 1975). Be- 
cause of the rapid changeover from a convec- 
tive to a radiative regime, our conclusions re- 
garding convection may be sensitive to the de- 
tails of the photospheric models, and at present, 
these are somewhat uncertain. Thus, although 


convection appears to be unimportant, this is 
not a robust conclusion, particularly if inhomo- 
geneities are present. 

Several authors have found it convenient to 
distinguish convective theories as being “local” 
or “nonlocal.” The local theories are con- 
structed so that, at a given position in a star, 
all quantities depend on local values of other 
variables, while nonlocal theories depend on 
values of variables at distant points in the at- 
mosphere. By all odds, the most widely used 
treatment of convection is the local mixing- 
length theory (Bohm-Vitense, 1958). This the- 
ory has the advantage of ease of computation, 
a free parameter which can, if necessary, be ad- 
justed to match observations, and comparabil- 
ity with past research. In the mixing-length the- 
ory (cf. Mihalas, 1978), the convective flux is 

F c = QC p vT(v- v E ) l/H , (7-5) 

where v is the average velocity of a bubble of 
gas, V = d In T/d In P in the ambient atmo- 
sphere, V £ refers to the value of V in the 
convective element, £ is the mixing length (the 
average distance traveled by a bubble before it 
loses all its energy), and H is the pressure scale 
height. The average velocity v is found by set- 
ting the kinetic energy of the bubble to one-half 
the work done by the buoyant force. The other 
half goes to overcome friction by other bub- 
bles. The gradients V and V E are calculated 
from considerations of the energy loss by the 
bubble during its lifetime and depend on as- 
sumptions regarding optical thickness of the 
bubble and its rate of radiation. To the present 
time, all significant sets of atmospheres for red- 
giant stars have been based on radiative equi- 
librium, or if convection has been included, it 
has been handled by the local mixing-length 
theory. Improvements made in or suggested for 
the standard mixing-length theory include tak- 
ing account of the optical depth of the convec- 
tive bubble (Henyey et al., 1965) and averag- 
ing the opacity across horizontal layers, rather 
than averaging the temperature and then calcu- 
lating the opacity (Deupree, 1979). 
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The question of the mixing length, or the ra- 
tio of mixing length to pressure scale height 
( £ / H ), has provoked considerable study. Most 
often one simply sets l /H = 1.0 or 1.5. For 
cool giant stars, surface convection is of such 
minor importance that £ ///apparently cannot 
be fixed from observations. (Indeed, it is not 
even certain whether convection can be treated 
by a local theory.) It is possible that the shape 
of the flux peak at 1.67 jtm (Bell et al., 1976b) 
or the profile of a strong LTE line such as Ca 
I 6572 or CO (Carbon et al., 1976) might be 
a valuable diagnostic feature for temperature 
structure and hence convection. In a different 
approach toward determining the scale length 
of convection, B5hm-Vitense and Nelson (1976) 
note the near coincidence in the HR diagram 
of: (1) the boundary line separating convective 
and nonconvective stars, and (2) the red edge 
of the Cepheid instability strip and use this fact 
to fix the mixing length at H, the pressure scale 
height. A later study (Bohm-Vitense and 
Dettman, 1980) demonstrates that chromo- 
spheric emission is observed by IUE only for 
giants cooler than F2, and by implication, this 
denotes the onset of convection. A comparison 
of IUE fluxes with model predictions (Bohm- 
Vitense, 1982) appears to confirm this conclu- 
sion. 

Based on scaling the quantity (V - V a J, 
which is related to the depth of the convective 
zone, Schwarzschild (1975) speculated that 
enormous granules (~ 10 8 km) might exist on 
Betelgeuse. Although Bohm-Vitense and Nel- 
son (1976) found no evidence for such giant ed- 
dies, such features may not be entirely theoret- 
ical. High-resolution spectrophotometry of the 
Ca II K line in a Boo (Chiu et al., 1977) ap- 
pears to require two states of excitation of the 
chromosphere. The normal state gives a mass 
loss dM/dt<lO~ 9 Mq /yr, whereas the excited 
states gives dM/dt ~ 8 X 10' 9 A/©/yr. The 
authors wonder whether the excited state arises 
from the appearance at the surface of a giant 
convective element. The argument for a two- 
component atmosphere is strengthened by a 
study of CO, which shows that no single- 
component atmosphere can match observations 


of both the wings of the Ca II infrared triplet 
(Ayres and Linsky, 1975) and the CO lines 
(Heasley et al., 1978). Of greater relevance here 
is Betelgeuse, whose Ca II and Mg II lines ap- 
pear to require an inhomogeneous atmosphere 
(Basri et al., 1981). Finally, observations of a 
regular variation in the polarization of this ob- 
ject (Hayes, 1984) strongly suggest rotation of 
a “hot spot” (Schwarz and Clarke, 1984; 
Clarke and Schwarz, 1984) — perhaps a large 
convective cell or a magnetic complex. 

Nonlocal theories of convection have not yet 
been widely adopted for cool star atmospheres. 
It may be, however, that the application of the 
most successful of these theories will constitute 
a major field of progress during the coming 
decade. A pioneering, nonlocal theory of con- 
vection (Spiegel, 1963) made an analogy with 
the equation of radiative transfer, the “source 
function” being chosen so that, in a uniform 
medium, this “generalized mixing-length” 
theory gave the same flux as the standard 
theory. In a “diffusion approximation”, this 
theory was later applied, with gray opacities, 
to the Sun (Travis and Matsushima, 1973a) and 
to G and K dwarf stars (Travis and Mat- 
sushima, 1973b). For his study of yellow giants, 
Parsons (1969) introduced a nonlocal element 
into the standard mixing-length theory by 
averaging the convective velocity over the en- 
tire trail of the convective elements. Convec- 
tive thermal plumes in the terrestrial at- 
mosphere provided an analogy for a totally dif- 
ferent nonlocal theory (Ulrich, 1970a), which 
was then applied to the Sun, where significant 
overshooting into optically thin regions was 
found (Ulrich, 1970b). A critical review of these 
theories and numerical comparison of their 
predictions for the Sun was made by Nordlund 
(1974). With a correction in the Ulrich theory, 
its similarity to the schemes of Parsons was 
demonstrated. The low value of £ /H required 
by the Travis-Matsushima formulation to 
match solar limb darkening was explained as 
a result of excessive overshoot. By a proper 
choice of such free parameters as the mixing 
length, the methods of Parsons, (modified) 
Ulrich, or even the simple mixing-length theory 
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can be made to yield similar results and may 
confidently be used to calculate models of the 
optically thin photospheres of main-sequence 
stars and giants down to KO. In deeper layers, 
in which convection is more important, the 
theories offer qualitatively correct but 
parameterized descriptions of convection. 

Nonlocal theories, particularly those of 
Spiegel (1963) and Ulrich (1970a, 1970c), were 
compared by Ulrich (1976), who showed that 
these differ principally in the way in which the 
averaging is done. The length over which 
averages are taken — the diffusion length — thus 
becomes a significant parameter. The choice of 
a particular kernal function for the integration 
leads to the more tractable differential equa- 
tion found by Travis and Matsushima (1973a, 
1973b). 

Another approach is the anelastic convec- 
tive theory (Latour et al., 1976), which borrows 
from meteorology the idea of greatly restricting 
the allowed acoustic modes. Horizontal modes 
are also simplified to permit a more detailed ex- 
amination of the large-scale vertical structure. 
In subsequent papers, the anelastic convective 
theory has been applied to convection in an A 
star, for which the convective zones are well 
separated and convection is relatively modest 
(Toomre et al., 1976; Latour et al., 1981). 

All of the foregoing theories, whether local 
or nonlocal, are based on a “one-stream” 
model of convection in which convective ele- 
ments (bubbles or plumes) move relative to a 
static ambient atmosphere. A more powerful, 
but much more complex, scheme would be a 
“two-stream” model. A number of such semi- 
empirical models have been constructed to in- 
terpret velocity and brightness inhomogeneities 
in the solar photosphere. A more or less deduc- 
tive two-stream model of convection was de- 
rived by Ulrich (1970a, 1970c), and a fully 
deductive model was derived by Nordlund 
(1976). This latter method has been applied to 
convection in the Sun and cooler stars with ex- 
cellent results (Nordlund, 1980, 1982). It is 
possible to simulate realistically photospheric 
granular convection and to predict reasonable 
values of line shifts and widths, which in turn 


give rise to what is called macroturbulence and 
microturbulence. Overshooting of convective 
elements into superficial radiative layers is 
specifically allowed, and the resulting large 
temperature differences ( — 10 to 10 3 K) bet- 
ween the two streams in visible layers of the 
photosphere constitute one of the most 
dramatic differences between these results and 
those of the local mixing length theory. At large 
depths (r > 50), the predictions of the two- 
component theory and those of the local 
mixing-length theory are quite similar. The in- 
terplay of various terms in the energy balance 
between the hot (ascending) components and 
the cool (descending) components at various 
heights is complex and instructive. In the outer- 
most layers, where the radiation begins to 
escape directly to the surface, these temperature 
fluctuations are due principally to radiative 
terms — radiative heating from below, radiative 
exchange between the two components, and 
radiative losses to the surface — which implies 
that the temperature difference in the optical- 
ly thin layers depends on the horizontal scale 
of the convecting elements and emphasizes the 
crucial role of that quantity (Schwarzschild, 
1975, Bohm-Vitense and Nelson, 1976). Clear- 
ly, there is a hint of far-reaching potential in 
this two-stream theory, and we eagerly look 
forward to future advances (Nordlund, 1982, 
1985) and applications to red-giant stars. 

A fascinating study of the effect of certain 
improvements in the mixing-length theory and 
a comparison with other theories has now been 
made for the Sun (Lester et al., 1982). Two im- 
provements were incorporated into the 
ATLAS6 treatment of convection: (1) the 
proper horizontally averaged opacity (Deupree, 
1979) was taken, and (2) a depth-variable mix- 
ing length was employed. The results are sum- 
marized by the authors as follows, in compar- 
ing the new calculations with otherwise identical 
calculations for the Sun without these modifi- 
cations. 

1. Convection (in the improved treatment) 
transports a smaller fraction of the total 
flux. 
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2. The convection zone is narrower, and 
F /F, is smoother. 

conv tot 

3. The temperature is higher in the convec- 
tive region. 

4. The model is brighter in the region 1500 
to 2000 A. 

5. Small differences arise in the visual spec- 
trum. These differences result almost en- 
tirely from the use of the horizontally av- 
eraged opacity. Compared to other the- 
ories, there are also differences, and 
these result from both the horizontally 
averaged opacity and the variable mix- 
ing length (Deupree and Varner, 1980). 

6. Compared to nonlocal mixing-length the- 
ories, the new calculation has less con- 
vection in the upper photosphere, and 
this leads to these models being hotter be- 
low t r = 0.6. 

7. Compared to Nordlund’s (1976) two- 
stream model, the new results agree well 
above r = 15, but are much hotter below 
that depth. 

8. Compared to empirical solar models, 
there is good agreement above log r - 

0.4, but the new results are too warm be- 
low that point. 

The Sun is a natural testing ground for new 
theories, and it may be some time before better 
theories are applied to the red-giant stars which 
form the subject of interest here. At present, 
the entire field of convection is fragmented and 
in turmoil (as if convecting?), making a coher- 
ent review next to impossible (cf. Marcus et al., 
1983). Worse yet, for our purposes, application 
of the available theoretical work is almost ex- 
clusively to the Sun, for which high-quality ob- 
servations of limb darkening, granulation, and 
line profiles can be obtained and an empirical 
model is available. Improvements in convective 


theory for red-giant stars will likely be bor- 
rowed from applications in warmer stars. Two 
comments are relevant here: (1) inhomogenei- 
ties are likely to be relatively much more im- 
portant in red-giant stars than in the warmer 
more compact stars in which convection is usu- 
ally studied; and (2) because red giants are 
variable and thermal models are crude, the 
power of this combination to constrain or guide 
the theory of convection is weak. 

In summary, we note that convection and 
inhomogeneities might play important roles in 
at least the following ways. 

1 . Because concentrations of certain mole- 
cules (H 2 0, H 2 , HCN, TiO, and CaCl) 
are extremely temperature-sensitive, their 
concentrations must increase over cooler 
regions of the photosphere. One can 
imagine a feedback mechanism through 
which increased molecular concentra- 
tions further cool the atmosphere 
(Gustafsson, 1981). 

2 . The subphotospheric convective zone 
presumably supplies the mechanical en- 
ergy to produce a chromospheric temper- 
ature rise and perhaps promote mass loss 
(cf. Schmitz and Ulmschneider, 1981; 
and de la Reza and Goldberg, both this 
volume). 

3 . Grain formation might occur, or occur 

more readily, in or over cooler photo- 
spheric regions (Lucy, 1976; Deguchi, 
1980; Draine, 1981, 1984; Woodrow and 
Auman, 1982; see also Lefevre, this 
volume). 

4 . Microvelocity fields deduced from high- 

resolution spectra — or their interpreta- 
tion — would be affected (Dravins et al., 
1981; Nordlund, 1982). 

5 . Departures from LTE might be in- 
creased, perhaps dramatically, in an in- 
homogeneous atmosphere. 
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6 . The conversion of theoretical quantities 
(such as T eff ) to observed quantities 
(such as colors) might be affected, as has 
already been noticed for F stars (Nelson, 
1980). 

INTERLUDE 

How well do the predictions of the theoreti- 
cal models of the previous section compare with 
the observations of real stars? A moment’s re- 
flection suffices to remind us, of course, that 
these thermal models will yield only a back- 
ground continuous spectra and symmetric ab- 
sorption lines, and there is no possibility that 
thermal models can reproduce the dazzling ar- 
ray of observations presently available (as de- 
scribed, for example, by the Quercis in Chap- 
ters 1 and 2). Much more work is necessary 
before we can model, in a self-consistent 
fashion, such nonstatic and nonthermal 
phenomena as turbulence, departures from 
LTE, pulsation, and mass loss (see Figure 7-1). 
Progress toward these goals is reported in later 
chapters. It is only fair, however, to give the 
thermal models a chance to demonstrate 
whether they can explain even what they were 
meant to explain. In addition, we seek spectral 
features which are particularly sensitive to a 
single property of the atmosphere and may 
therefore help to diagnose the physical and 
chemical properties of the star. At the outset, 
we set forth our most difficult problem: how 
does one test the goodness of a comparison bet- 
ween theory and observation? As will become 
obvious from several examples, we will stum- 
ble over this problem again and again. 

SPECTRAL DISTRIBUTION 
OF ENERGY FLUXES 
(Photometry, Spectro- 
photometry, and Colors) 

Oxygen-Rich Stars (M and S Stars) 

Although the models of Gustafsson et al. 
(1975) and Bell et al. (1976a) are intended for 
G and K giant stars, we mention them here be- 


cause of their obvious implications for cooler 
stars and the paucity of such comparisons for 
M, S, and C stars. The excellent overall match 
of observed spectrophotometry and narrow- 
band photometry with predicted fluxes from 
these models confirms that the line-blanketing 
and thermal structure of present models of G-K 
giants is essentially correct (Gustafsson and 
Bell, 1979). This impressive work has reached 
a new level of sophistication for stars of inter- 
mediate temperature. Unfortunately, the red- 
giant stars of interest in this volume constitute 
a rather different class of stars. 

A persistent flux excess shortward of 4000 
A was carefully scrutinized by Gustafsson et al. 
(1975), who conclude that, although massive 
chromsopheres, temperature inhomogeneities, 
and departures from LTE cannot be ruled out 
as sources for the discrepancy, the most likely 
interpretation is neglected ultraviolet 
opacities — from either metal lines, molecular 
lines, or some continuous absorber. A violet 
discrepancy is not unique to these models, but 
has been noted in other models as well. For ex- 
ample, Peytremann (1974) found that his 
models produced too much flux in the violet, 
a feature he also suggested was due to insuffi- 
cient lines there. In her models of F dwarf stars, 
Bohm-Vitense (1982) finds the same cir- 
cumstance and suggests that it arises from either 
an unaccounted-for continuous opacity source 
or a substantially increased convective flux. A 
similar problem appears in the M giant models 
of Johnson et al. (1980). (See Steiman-Cameron 
and Johnson, 1985.) To investigate the possi- 
ble luminosity dependence of the relative flux 
peak at 1 .67 /an, for which there is some obser- 
vational evidence, ODF spectra of CN and CO 
were computed and found to be in fair agree- 
ment with the observations (Bell et al., 1976b). 
These conclusions are only tentative until ad- 
ditional observations at many frequencies, as 
well as improved theoretical spectra including 
the contributions of TiO and H 2 0, are 
available. Limb-darkening coefficients for the 
set of models have also been calculated (Man- 
duca et al., 1977) for the wavelengths of the 
UBV and uvby bandpasses. 
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CALCULATED/OBSERVED FLUX 


Based on Gustafsson et al. (1975) and Bell 
et al. (1976a) atmospheres with 4000 < T ff < 
6000 K, Bell and Gustafsson (1980) calculate 
the visual surface brightness (VSB) and study 
theoretically its dependence on color, gravity, 
and metallicity. The total effects of all these 
factors are smallest for the V-R color, and the 
theoretical relation between VSB and V-R 
agrees well with the empirical relation from 
stars with known angular diameters (Barnes 
and Evans, 1976; Barnes et al., 1978; Eaton and 
Poe, 1984). This agreement provides additional 
evidence that the background continuous spec- 
trum, including the effects of line blanketing, 
is well predicted by these theoretical models for 
G and K giants. 

Spectrophotometry provides a somewhat 
stiffer test of models than does photometry, 
and the comparison of predictions from the 
Gustafsson et al. (1975) and Bell et al. (1976a) 
models with spectrophotometry of G, K, and 
M stars (Strecker et al., 1979) is of considerable 
interest. Water vapor, but not TiO, was 
included in the synthetic spectrum, although 
neither was included in the calculations of the 
models. Figure 7-12 shows the comparison for 
ft And (M0 III) and (3 Peg (M2 II-III) (Man- 
duca et al., 1981); the agreement over most of 



Figure 7-12. Comparison of observed spectro- 
photometry of early M giant stars & And and 
(3 Peg (Strecker et al., 1979) and predictions 
from the Gustafsson et al. (1975) and Bell et 
al. (1976a) models identified in the figure (from 
Manduca et al., 1981). 


the range from 1.2 to 5.5 fim is better than 10 
percent. (Fits for warmer models are generally 
even better.) The largest discrepancy in ft And 
occurs in the first overtone bands of CO at 2.5 
pa n. While that particular problem could in 
principle be remedied by altering the CNO 
abundance, the ratio of the fundamental and 
the first overtone of CO cannot be so easily cor- 
rected. The fit for ft Peg is slightly worse, but 
this comparison should not be pushed too far, 
since the authors suspect that the star may be 
cooler than any of their models. The authors 
also note that effective temperatures deduced 
from model fitting are in satisfactory agreement 
with other temperature estimates, and angular 
diameters inferred from fitting at 3.5 fim are 
in good agreement with measured values of 
those deduced from the VSB relation and V-R 
colors. 

Comparison of computed fluxes from a grid 
of SM models (Johnson, 1974) with infrared 
spectrophotometry of 24 cool giant stars (K0 
through M9) taken from above the Earth’s at- 
mosphere (Scargle and Strecker, 1979) furnishes 
some food for thought. As a matter of fact, the 
predicted fluxes often agree rather well with the 
observations. Yet some of this agreement must 
surely reflect the insensitivity of the fluxes to 
the precise structure of the models because SM 
models have outer photospheres that are 
generally too warm. More surprisingly, many 
of these “models” were interpolated or ex- 
trapolated where models did not exist in the 
original grid. Furthermore, there is a strange 
clustering of stars with T eff - 3500 K, which 
appears to be an artifact of the available set of 
models, and the temperatures from several of 
the models are lower than those found for some 
of the same stars by other methods (Gustafs- 
son, 1981). All of this is a bit disquieting 
because it forces us to admit that, unless the 
effective temperature is constrained rather care- 
fully from other evidence, observations, 
especially for variable stars, can be matched 
reasonably well by crude models. Better obser- 
vations , especially calibrated spectrophotome- 
try of both the visual and infrared regions taken 
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at the same time, are necessary to test even ther- 
mal models . 

Tsuji (1976b) found good agreement be- 
tween infrared fluxes from his models of M su- 
pergiants and accumulated photometry and 
spectrophotometry of Betelgeuse. The fit was 
improved if the composition was taken to be 
C-poor and N-rich, rather than solar. (That the 
atmosphere of Betelgeuse is actually C-poor 
and N-rich has recently been determined 
through a careful study using spectral synthesis 
(Lambert et al., 1984).) A relatively high micro- 
turbulent velocity of 6 km/s was used; this was 
not taken as evidence for supersonic turbulence, 
but was attributed to deficiencies in the 
VAEBM method for molecular line opacities. 
The best fit to the infrared fluxes (1 to 5 pan) 
gave an effective temperature of 3900 ± 150 
K for Betelgeuse (M2 lab). The implied angular 
diameter of 0.40 milliarc-sec is, however, lower 
than most other determinations. No resolution 
of this discrepancy was offered, but Tsuji noted 
that this higher temperature resolves a problem 
with the K I line (Goldberg et al., 1975; see also 
Goldberg, 1984). (Could it be that most earlier 
observations detected the larger radii in TiO 
bands?) Subsequently, Tsuji (1978b) demon- 
strated that scattered light from an optically 
thin dust shell had negligible effect on the col- 
ors , but it might appreciably increase the mea- 
sured angular diameter , which might resolve the 
Betelgeuse discrepancy. The strength of the TiO 
bands (Collins, 1974) implied by the variation 
of angular diameter of Betelgeuse across them 
appears to require a lower temperature, how- 
ever, than 3900 K (Balega et al., 1982). The cu- 
rious behavior of the angular diameter of 
Betelgeuse is quite complex and is not yet un- 
derstood (Goldberg, this volume). 

In a pioneering extension of model atmo- 
spheres to cooler stars, Tsuji (1978a) saw suffi- 
cient agreement between fluxes from his models 
and available photometry and spectrophotom- 
etry of K and M giant stars to define a new scale 
of effective temperature for these objects. A 
careful scrutiny of these comparisons reveals 
that, in the infrared, the slopes of the flux 
curves of all models are quite similiar, and most 


of the temperature sensitivity comes from wave- 
lengths less than 1 pa n. Here, however, none of 
the predicted fluxes fit well because of distor- 
tions by molecular bands. The fitting is there- 
fore fairly subjective. Nevertheless, if R-I col- 
ors are used as a baseline, there is reasonable 
accord between TsujFs scale of effective tem- 
perature and the scale deduced from the angu- 
lar diameters then available. Finally, Tsuji 
pointed out that colors and band strengths 
could well be understood on the basis of this 
new temperature scale, but no information 
about variations in abundances among the 
CNO elements could be inferred. 

As additional values of stellar angular di- 
ameters became available, especially from ob- 
servations of lunar occultations, it became pos- 
sible to define an empirical temperature scale 
for red giants from K0 to M6 (Ridgway et al., 
1980b). Although the scatter among individual 
stars is rather high and certain large discrepan- 
cies exist even for stars for which the formal 
errors are small, this calibration of effective 
temperature against spectral class or Wing color 
temperature ( T c ) has become the criterion 
against which all theoretical results are tested. 
From time to time, more values of angular di- 
ameters and temperatures have been obtained, 
but no basic alteration of the Ridgway et al. 
scale has been necessary. 

An additional link between models and 
observations was forged by Tsuji (1981a), who 
used his atmospheres and the method of in- 
frared photometry (IR) (Blackwell and Shallis, 
1977) to define a new independent temperature 
scale for M giant stars. Basically, the IR 
method is a clever application of the relation 
between observed flux (/ x ), flux emitted at the 
stellar surface (F x ), and the stellar angular 
diameter, in which the monochromatic flux 
from an appropriate model is substituted for 
the emitted stellar monochromatic flux (F x ) 
and is used with the observed monochromatic 
flux (/ x ) to obtain the angular diameter. With 
this value of the angular diameter and the in- 
tegrated observed flux (/), the integrated emit- 
ted stellar flux (F) and hence the effective 
temperature can be derived. The derived effec- 
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tive temperature depends only rather weakly on 
the model used. Furthermore, since the ratio 
of IR flux to total flux varies roughly as r ff 3 , 
the method is quite powerful in determining the 
temperature even in the presence of imprecision 
in the models from missing opacities, inexact 
opacity treatment, or imperfect match in chemi- 
cal composition. The necessary equations have 
been combined into an especially convenient 
form (Blackwell et al., 1980). From photomet- 
ric measurements on 62 giant stars from spec- 
tral class K1 through M6, Tsuji (1981a) used 
L-band photometry for the IR fluxes to deduce 
effective temperatures. The resultant tempera- 
ture scale agrees quite satisfactorily with that 
deduced from angular diameters (Ridgway et 
al., 1980b). A minor difference is that Tsuji 
found evidence for a plateau at F ff = 3800 K 
for stars of spectral types MO through M2 and 
speculated that, in this spectral range, luminosi- 
ty might affect the spectral classification in such 
a way that stars of identical temperature but 
differing luminosities were assigned different 
spectral classes. (The reader may detect a minor 
inconsistency here, because the values of F eff 
calculated from the models do not show a 
gravity dependence.) Although a slowing of the 
steady decrease of T ff with spectral type can- 
not presently be excluded by angular diameter 
data, Keenan (1982) has demonstrated that 
there is a steady progression of color and spec- 
tra of Fe, Cr, and TiO lines through these spec- 
tral classes. 

Another crucial issue, noted but not resolv- 
ed by Tsuji (1981c), arises from the masses of 
the M giants inferred from various luminosity 
measurements and his effective temperatures. 
In the resultant HR diagram (Figure 12, Tsuji, 
1981c) the K-M giants scatter along the theo- 
retical evolutionary track for a star of 3 solar 
masses, yet present knowledge of the luminosity 
and initial mass function appears to establish 
that most of these stars must have masses near 
1 solar mass (Scalo et al., 1978; Scalo and Mil- 
ler, 1979). 

Based on the Johnson et al. (1980) models 
of K and M giants described earlier, Piccirillo 
et al. (1981) calculated Wing color temperatures 


and compared the relation between T c and F ff 
with that obtained from previous work. The 
original calibration of T q against F ff is semi- 
empirical (Ridgway et al., 1980b). As is clear 
from Figure 7-13, in their respective ranges of 
overlap, the models of Gustafsson et al. (1975) 
and Bell et al. (1976a), of Tsuji (1978a), and 
of the Indiana group (Johnson et al., 1980) 
agree well with the semiempirical relation. This 
agreement of the three sets of theoretical 
models with the semiempirical relation con- 
stitutes the most convincing evidence that 
theoretical models of early M giants are basical- 
ly correct . Piccirillo et al. (1981) further 
demonstrate that the difference between T q 
and T ce is due to a multitude of weak TiO 
lines, unseen even at fairly high spectral resolu- 
tion, which depress the stellar flux by signifi- 
cant amounts in the region of the short- 
wavelength continuum filter at 7540 A (Wing, 
1981). 



Figure 7-13. The relation between Wing color 
temperature and effective temperature for K 
and M giants. Solid line is the semiempirical re- 
lation (Ridgway et al., 1980b); symbols repre- 
sent results obtained from model calculations 
by three different groups (Gustafsson et al., 
1975, and Bell et al., 1976a; Tsuji, 1978a; John- 
son et al., 1980). Note that the points calculated 
without TiO absorption appear to follow the 
equality line. References are given in the text 
(from Piccirillo et al., 1981). 
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Results from the first comprehensive com- 
parison of colors from models of K-M giant 
stars and observations are also quite encourag- 
ing (Steiman-Cameron and Johnson, 1985). 
Magnitudes on the VRIJHKL broadband 
system and on the wing narrowband system 
have been computed for the models of Johnson 
et al. (1980) and compared to observations. 
Figure 7-14 shows a comparison of effective 
temperature against the TiO band strength in- 
dex computed from the Johnson et al. models, 
compared with the relations of Ridgway et al. 
(1980b) and Tsuji (1981a). We note that the 
Ridgway et al. curve is semiempirical (the best 
fit by eye through the observed points shown), 
the Tsuji curve is also semiempirical since it uses 
T „ deduced by the IR method and TiO 
measured for the corresponding spectral 
classes, and the Johnson et al. curve is purely 
theoretical. The agreement is quite satisfactory, 
certainly within the large observational scatter, 
down to a temperature of 3500 K, below which 
the theoretical TiO band strengths are too 
weak. Steiman-Cameron and Johnson (1985) 
have considered several possible effects to ex- 
plain the discrepancy at lower temperatures, in- 
cluding dissociation of TiO by chromospheric 
heating, increasing effects of sphericity as the 
luminosity increases and the surface gravity 
decreases toward cooler stars, changed thermal 
structure due to possible carbon enhancement 
in the cooler stars, and a variation of radius 
with wavelength. (A chromospheric effect ap- 
parently exists (Steiman-Cameron et al., 1985).) 
In any case, the excellent accord for early M 
stars is reassuring. 

Carbon Stars and S Stars 

What of the fascinating peculiar red-giant 
stars — types N, S, and R? Most model calcula- 
tions are so recent that time has not allowed 
careful comparison with observations, and such 
comparisons remain a pressing need. The as- 
surance we could give earlier that predictions 
from classic thermal models for early M giants 
were in reasonable accord with observations of 
the continuous flux, and that reasonable values 


of angular diameters and temperatures could 
be deduced on the basis of model fitting, can- 
not be repeated for the cooler stars. 

The French group (Querci and Querci, 1976) 
compare two of their triple-alpha-enriched 
carbon-star models with the observed broad- 
band colors of 19 Psc and UU Aur. The fit for 
a model with T , = 3600 K to UU Aur is fair, 

eft 

and the theoretical predictions fit better than 
the blackbody curve, especially in the 0.8- to 
1.6-/un region. The flux curve for 3000 K 
matches almost as well, but is disfavored by 
Querci and Querci because of a poor fit to the 
observed spectrum. All models predict excess 
fluxes in the blue part of the spectrum. 
Altogether, the comparison leaves much to be 
desired because the fit is only fair, and the tem- 
perature deduced from the best fit — T ff = 
3600 K for UU Aur for the triple-alpha compo- 
sition — appears to be considerably too high 
(Ridgway et al., 1980a; Tsuji, 1981b). A lower 
temperature is obtained for a model based on 
CNO processing, but synthetic spectra are then 
in conflict with observations. One is hard 
pressed, as were Querci and Querci (1976), to 
select from a large number of choices the spe- 
cific causes for the relatively poor fits — possibly 
including at least incorrect chemical composi- 
tion, neglected opacities, and neglected convec- 
tion. 

In a bold extension of his earlier work on 
M giants, Tsuji (1981b) used the available mod- 
els of carbon stars (Querci et al., 1974; Querci 
and Querci, 1975a) and the L-band observa- 
tions of Noguchi et al. (1977) to deduce a 
temperature scale for N-type carbon stars by 
the IR photometry method described earlier. In 
this case, it was necessary to correct theoretical- 
ly the L-band observations for molecular ab- 
sorption due to HCN and CO. The resulting 
temperatures are in satisfactory agreement with 
the six obtained by lunar occultation angular 
diameters (Ridgway et al., 1980a). Tsuji (1981c) 
further demonstrates that several broadband 
colors of carbon stars can be understood on the 
basis of this effective temperature scale. 

Recently, NASA/Ames Research Center 
spectrophotometry of 10 non-Mira S stars and 
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Figure 7-14. The relation between effective temperature and spectral type (TiO band strengths) 
for K and M giant stars. Shown are the observations , the semiempirical curve which best 
fits the observations (Ridgway et al., 1980b), the semiempirical relation from IR photometry 
( Tsuji , 1981a), and the theoretical fit by the models of Johnson et al. (1980). (From Steiman- 
Cameron and Johnson, 1985.) 


fluxes from theoretical models (Johnson, 1982) 
have been compared in an attempt to calibrate 
a temperature scale for S stars (Augason et al., 
1985). Effective temperatures have been ob- 
tained both by model fitting and by the method 
of IR photometry; results from the two meth- 
ods are in good agreement. An example of the 
observed and calculated spectrophotometry for 
the warm S star, HD 35155, is shown in Figure 
7-15. The relative strengths of the fundamen- 
tal and the first overtone bands of CO do not 
agree well (a problem we have met before). The 
keen observer will also note the unusual flat- 


topped region near the maximum of the H” 
peak at 1.67 ^m, which differs strikingly in ap- 
pearance in S and C stars. If the observations 
were extended to shorter wavelengths, an ex- 
cess of flux in the visual and ultraviolet spec- 
tral regions from the model would also be no- 
ticed, and comparisons of the theoretical flux 
curves to the observations of cooler S stars are 
not as good as the fit shown. Also, ZrO is not 
included in these models, and ZrO bands are 
prominent in the comparison. Although judg- 
ment is left to the reader, it is clear that prog- 
ress is being made in understanding even these 
complex objects by thermal models. 
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Figure 7-15. Comparison between spectropho- 
tometry of an early S star and predicted flux 
from an Indiana model (from Augason et al. , 
1985). 


A similar study has been made of the car- 
bon star, Y CVn (C5,5), based on extensive 
photometry (Goebel et al., 1980). A good fit 
to the entire spectral range 0.5 to 5.0 /im is ob- 
tained for theoretical fluxes from models 
(Querci et al., 1974; Querci and Querci, 1975a). 
Not only is the continuous spectral energy dis- 
tribution well matched, but the molecular band 
strengths of CN and C 2 are satisfactory. 
Whether the effective temperature is correct, 
is not yet clear, however, and one must be cau- 
tious of extrapolating general conclusions from 
an isolated case. 

Broadband photometry, which is often the 
only data available on a given star, carries a 
minimum amount of information, and the 
power of observations of one star to constrain 
the parameters of models is very weak. The 
power of photometric observations of many 
stars with differing properties to constrain a 
series of models is considerably enhanced. Nar- 
rowband photometry and spectrophotometry 
are considerably more powerful tools. Al- 
though we await eagerly such complete spectro- 
photometry for a wide variety of M, S, and C 
stars, we must be cautious, because unlike 
warmer stars, essentially all stars of interest 
here are variable. Combining observations 
taken at random phases — sometimes of neces- 


sity — introduces a significant uncertainty into 
the results. 

We conclude that thermal models of the 
classic variety are able to reproduce the con- 
tinuous flux curves of giant stars of solar 
composition in spectral classes MO through M5. 
Progress is also apparent for S stars, especially 
the warmer ones. Testing of C-star models is 
just beginning. Many further improvements on 
models are necessary. Continually updated 
opacities are needed. More molecular lines may 
be needed for familiar molecules (the huge 
numbers of weak lines), and opacities from ad- 
ditional molecules should be included. For cool 
M stars, the lack of an accurate opacity for 
H 2 0 is perhaps the single greatest impediment 
to further progress, but VO is also important. 
In S stars, ZrO may be important. In cool car- 
bon stars, one must add HCN (Eriksson et al., 
1984; Tsuji, 1984; J0rgensen, 1985), C 2 H 2 , 
and perhaps C^H, C y and SiC 2 . Perhaps we 
should look at SiO, CS, SiS, FeH, HC1, A1C1, 
and CaCl — the list could be long. 

It is also in these coolest M, S, and N stars 
(especially the M stars, according to Scholz and 
Tsuji (1984)) that the effect of sphericity is ex- 
pected to become important as the gravity de- 
creases with the increase in radius during these 
late stages of stellar evolution. It is encouraging 
to note the progress being made by the Heidel- 
berg group and the Japanese workers. Some of 
our fitting or lack of fitting could be pro- 
foundly changed as a result of sphericity. How- 
ever, a gulf still exists here. The model makers 
who use plane-parallel models have the exten- 
sive opacity data and, understandably, continue 
their zealous research with these models. Surely 
there is much room for that. The spherical 
models have generally been converged with rel- 
atively crude opacities. It is time for a merging 
of forces — of those that have the good opacities 
with those that have the good spherical codes. 

SPECTRAL LINES 

There are at least three vital questions, the 
answers to which will be essential as we con- 
sider to what extent line profiles from classic 
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thermal models match observations of real 
stars. How can one judge quantitatively good- 
ness of fit? Can one isolate the specific reason 
for the lack of perfect fit? What can one deduce 
regarding physical conditions in the stellar at- 
mosphere from various lines? The questions are 
especially pertinent to red-giant stars, in which 
variability is known and departures from plane- 
parallel geometry, LTE, and homogeneity are 
expected. Present progress in answering these 
questions for K giants can be learned from the 
detailed review of the conference on Arcturus 
(Trimble and Bell, 1981). Much less work has 
been done on the red-giant stars. In addition, 
spectral lines in supergiants are remarkably 
broad, and considerable macroturbulence must 
often be invoked to match observations. 

An example of attempts to deduce 12 c/ 13 c 
ratios in carbon stars will illustrate the problem 
in answering definitively the foregoing ques- 
tions. Values obtained for the well-observed N- 
type star, V460 Cyg, from excellent high-reso- 
lution spectra by several methods make an en- 
grossing story: (1) 100, from a pseudo-curve- 
of-growth and the isointensity method (Fujita 
and Tsuji, 1977); (2) 9, from the Minnaert semi- 
empirical formula (Climenhaga et ah, 1977); (3) 
32, from a curve of growth (Dominy et al., 
1978); (4) 24, from spectral synthesis and a 
model atmosphere (Johnson et al., 1982); (5) 
22-29, from a careful reapplication of the Min- 
naert formula (Little-Marenin and Little, 1984). 
The precise values are of less importance in this 
discussion than the fact that crude methods and 
atmospheres gave results as good (or bad) as 
those for better atmospheres. We are evident- 
ly in the awkward position of admitting that 
crude model atmospheres can sometimes both 
produce good fits to fine observations and per- 
mit the inference of useful quantities about the 
star. Apparently, fitting observed spectra is a 
necessary, but not sufficient, condition for ac- 
cepting a thermal atmosphere! Clearly one 
needs tougher tests for atmospheres: broader 
spectral intervals, lines from several atoms, 
ions, and molecules and lines spanning a larger 
range in region of formation (hence tem- 
perature). Otherwise, in the absence of fairly 


strong constraints, goodness of fit to average 
lines is a very weak criterion for the identity of 
the model with a real star. 

A useful illustration of the problems of de- 
ducing stellar parameters from a comparison 
of observed and synthetic spectra is provided 
by the pioneering attempt to determine H, C, 
N, and O ratios in UU Aur and Y CVn (Querci 
and Querci, 1976). Such comparisons, the best 
of which is illustrated for UU Aur in Figure 
7-16, provide some reassurance that the models 
generally have a correct thermal structure, but 
uncertainty in other parameters do not allow 
a more detailed test. For example, the models 
shown give T eff ~ 3600 K for UU Aur, yet this 
is much higher than temperatures inferred for 
similar N-type stars (Ridgway et al., 1980a; 
Tsuji, 1981b). Compositions enriched in C 
(called PL in Figure 7-16) fit generally better 
than those with CNO-cycle compositions 
(called DE), but line and band /-values are not 
sufficiently well known to extract more precise 
information concerning H, C, N, and O 
abundances. Although the influence of gravity 
on line profiles is quite small, it increases 
slightly toward the lower temperatures. Its in- 
fluence varies somewhat from band to band; 
for example, CO and CN behave oppositely, 
and each behaves differently for PL and DE 
compositions. Because only one value of the 
turbulence, v = 5 km/s, was used, little infor- 
mation can be deduced from it. Reasonable 
overall fits to observed spectra of both stars in 
three spectral regions could be obtained, and 
the lack of better matches can also be attributed 
to imprecision of molecular data as to deficien- 
cies in the atmospheric models. We will find 
other examples of this sort of impasse. 

In their comprehensive study, Gustafsson 
and Bell (1979) compared LTE synthetic spec- 
tra from their models for G and K giants to 
observed spectra for several wavelength bands; 
the overall fit is fairly impressive. 

An attempt at a more objective approach to 
the comparison of synthetic and observed spec- 
tra was applied to five carbon stars (Olson and 
Richer, 1979) to infer the values of X(CNO), 
turbulent velocity, and 13 C/ 12 C by minimizing 
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the residuals between observed spectra and syn- 
thetic spectra based on earlier models (John- 
son, 1974). Such differences depend critically 
on accurate values for wavelengths and other 
molecular parameters, and perhaps the method 
was pushed a bit beyond the precision of the 
data. Yet such methods appear to hold consid- 
erable promise for the future as better observa- 
tional and molecular data become available. 
Our present interest centers on the goodness of 
the fit. The authors note that their synthetic 
spectra fit very well the high dispersion (13 
A/mm) observations of TX Psc, Z Psc, X Cnc, 
UU Aur, and Y CVn in the region 7800 to 8300 
A for reasonable values of the parameters. Un- 
fortunately, however, no absolute criterion ex- 
ists for judging goodness of fit. One of the 
parameters determined by fitting is the micro- 
velocity. The values derived by Olson and 


Richer (1979) are all subsonic ( Vs = 5 km/s) 
and significantly lower than values obtained by 
other workers. 

We await the results of a more extensive 
study of 12 C/ 13 C in approximately 30 N-type 
carbon stars from infrared Fourier transform 
spectrometer spectra and synthetic spectra from 
unpublished atmospheres (Lambert, 1984; 
Eriksson et aL, 1985; Gustafsson, 1985). 

One of the more significant results to come 
from spectroscopy of peculiar red-giant stars 
concerns the discovery and strength of the 
quadrupole lines of H 2 in the infrared (cf. 
Lambert et al., 1973; Hall and Ridgway, 1977; 
Goorvitch et al., 1980; Bragg et al., 1982; Jen- 
nings et al., 1984). The 1-0 S(l) line of H 2 at 
2.1 fir n — the line most often studied because of 
its strength and its location in a clear region of 



Figure 7-16. Spectra of the carbon star UU Aur and synthetic spectra from theoretical models with 
compositions representing triple-alpha-enriched material, by the French group (from Querci and Querci, 
1976). 
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the spectrum — has been detected and its 
strength followed throughout the cycle in the 
two S Miras, Chi Cyg and R And (Hall and 
Ridgway, 1977; Hinkle, private communica- 
tion, 1983). The same line has been detected in 
the carbon stars, V CrB, U Cyg, RZ Peg, and 
WZ Cas, but is not present in eight other re- 
type carbon stars, most of which are irregular 
or semiregular variables (Johnson et al., 1983). 
In a parallel investigation, Tsuji (1983) found 
no trace of this line in any M giant star, even 
though he studied stars as cool as RZ Ari (M6), 
WX Vir (M7), and RX Boo (M8). These results 
are even more unexpected because the line is 
predicted to be strong by published model at- 
mospheres (Querci et al., 1974; Johnson, 1982; 
Tsuji, 1978a) in all stars of this temperature 
range. Several possible causes for the nonap- 
pearance of the H 2 line in carbon stars have 
been suggested (Johnson et al., 1983): (1) the 
temperature scale for carbon stars is too low 
by 500 to 800 K; (2) the models are in- 
accurate — the most important deficiencies be- 
ing the assumption of a plane-parallel geometry 
and the neglect of the opacities of polyatomic 
molecules; (3) H 2 might be dissociated by 
chromospheric heating; (4) H 2 might be 
underabundant because of some nonequi- 
librium process; (5) non-Mira carbon stars 
might be hydrogen-deficient. Several of these 
considerations apply to the M giant star as well. 
At the moment, the source of the discrepancy 
is unclear. The claim has recently been put for- 
ward that the strength of the 0-0 S(l) line of 
H 2 is correctly predicted in N-type carbon 
stars by atmospheres with HCN opacities 
(Eriksson et al., 1984). The outer atmospheres 
of cool carbon stars are undoubtedly complex 
and variable thermodynamic regimes, and we 
suspect that a full understanding may entail a 
lot more work. 

All of this is not directly relevant for testing 
thermal atmospheres because only the presence 
or absence of the line is at stake, and that alone 
is not a strong criterion for the goodness of the 
atmosphere. The H 2 quadrupole lines are ex- 
cellent diagnostics for the temperature and den- 
sity of the upper photosphere, however, and 


careful observations of these lines at various 
phases in several stars are urgently needed. 

DIAGNOSTIC FEATURES 

Having discussed in some detail the (relative- 
ly rare) comparisons of model predictions with 
observations, we now consider attempts to in- 
fer from observations, perhaps with the aid of 
model atmospheres, specific characteristics of 
a star. In other words, how can one infer for 
a particular star the value of a specified 
property? 

An excellent general overview of the 
strengths and weaknesses of traditional 
methods of using model atmospheres to infer 
(or diagnose) values of effective temperatures 
and surface gravity for stars of all types is that 
of Pecker (1973). Curves of growth are fre- 
quently used, and many good discussions of 
these are available, including both general 
references (Aller, 1960; Gray, 1976; Mihalas, 
1978) and specific discussions of the use of 
curves of growth with atmospheres. Required 
reading for this section is the comprehensive 
and sobering review of progress in determining 
effective temperature, surface gravity, micro- 
turbulent velocity, chemical composition, and 
mass for Arcturus (Trimble and Bell, 1981). 
Uncertain as these quantities are for Arcturus, 
they are far less certain for red-giant stars, for 
which observations are much less complete or 
accurate. Furthermore, red-giant stars are vari- 
able. Alksne and Ikaunieks (1981) have sum- 
marized some of the problems for carbon stars. 

Effective Temperature 

Effective temperature determinations for 
red-giant stars will continue to be based, into 
the foreseeable future, on measured angular 
diameters and integrated flux curves. Angular 
diameters are currently best obtained for the 
red-giant stars from lunar occultation (cf. 
Ridgway et al., 1980b), and this will continue 
to be the case for some time. However, in- 
terferometers operating in the infrared and later 
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in the visual— perhaps from space — will even- 
tually be developed, leading to a new break- 
through in our knowledge of stellar atmo- 
spheres. Although the concept of effective 
temperature will lose some of its meaning for 
the coolest and most extended stars if their radii 
are found to depend on wavelength (Bonneau 
et al., 1982; Balega et ah, 1982), the effective 
temperature will continue to be used. Because 
uncertainty in the temperature can mask several 
other problems, including bad models, its deter- 
mination is a matter of prime importance. We 
commend those who have tirelessly pursued this 
research over the years and encourage them to 
continue their efforts, improving them when 
possible by observing at several wavelengths. 

The uncertainties of both the lunar - 
occultation method and the reduction of such 
data have been carefully described (Ridgway et 
al., 1980a). We reiterate here a warning for the 
observations of red-giant stars, most of which 
are variable . Broadband colors of the cool 
carbon star, AQ Sgr, measured as part of a ma- 
jor photometric program (Walker, 1979), for 
example, are noticably different from re- 
measurements of the same star at the time of 
a lunar occultation (Walker et al., 1979). For 
determining effective temperature, therefore, 
the photometry for these variable stars should 
be obtained at the same time as the angular 
diameter . 

Effective temperatures have been obtained 
by the method of infrared photometry (Black - 
well and Shallis, 1977; Blackwell et al., 1980) 
for both M giants (Tsuji, 1981a) and cool car- 
bon stars (Tsuji, 1981b). Values for T ff de- 
duced in this manner are in good agreement 
with those obtained from angular diameters as 
described earlier. As infrared spectrophotom- 
etry becomes available, temperatures can be ob- 
tained from an application of the IR method 
at several wavelengths, thus reducing any possi- 
ble error. Such a project should be given high 
priority since it can be done with existing equip- 
ment and models. 

A paper which calibrates the relation of the 
visual surface brightness (VSB) to (V-R) color 
and claims to extend this relation to stars of all 


spectral types and luminosity classes (Barnes et 
al., 1978) has already been noted (see also 
Eaton and Poe, 1984). Although primarily con- 
ceived to yield angular diameters, the method 
also gives effective temperatures when bolo- 
metric corrections are available. The validity of 
the method for the red-giant stars of this vol- 
ume has not been securely established because 
of the paucity of angular diameters and the 
variability of these stars. The VSB method has 
been used to determine effective temperatures 
for R stars (Dominy, 1984) and has been reex- 
amined for N-type carbon stars (Johnson, pri- 
vate communication, 1984). 

While model fitting is commonly used for 
deducing effective temperatures of hot stars, its 
application to red-giant stars has been delayed 
by the scarcity of both realistic models and in- 
frared spectrophotometry. A single published 
example employs the N-type carbon star, Y 
CVn, for which a temperature of T ff = 2600 
K was inferred by fitting fluxes from model at- 
mospheres (Querci et al., 1974) to the observed 
infrared flux curve (Goebel et al., 1980). Ef- 
fective temperatures for 10 non-Mira S stars— 
ranging from 2800 to 3600 K— have also been 
obtained from infrared spectrophotometry by 
both model fitting and IR photometry (Auga- 
son et al., 1985). As additional fluxes from 
(present and future) theoretical models are 
calculated and more extensive spectrophotom- 
etry becomes available, model fitting should 
prove to be a powerful tool for deducing ef- 
fective temperatures. 

Color temperatures , such as (R + I) - 
(J + K), calibrated by comparison with M giants 
(Mendoza and Johnson, 1965) or by black- 
bodies (Scalo, 1976; Walker, 1980), have been 
used for carbon stars, and certainly have their 
utility today in the absence of more fundamen- 
tal calibrations. Narrowband photometric col- 
ors (cf. Baumert, 1972) may yield more precise 
temperatures. Additional insights on the rela- 
tion of colors and effective temperatures for 
carbon stars is given by Alksne and Ikaumieks 
(1981), Tsuji (1981c), and Alksne et al. (1983). 
All of these authors find broadband colors to 
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be generally well correlated with effective tem- 
peratures deduced for cool carbon stars. An op- 
posing view, based on studies of carbon stars 
in our Galaxy and the Magellanic Clouds 
(Cohen et al., 1981), holds that infrared broad- 
band colors are more indicative of composition 
than of effective temperature. What is really 
needed, of course, is a calibration of either 
broadband or narrowband colors (Ridgway et 
al., 1980b; Wing, 1981) in terms of effective 
temperatures for all types of red-giant stars, so 
that these colors may be conveniently used to 
gage the effective temperature. 

A comprehensive scheme for determining 
the atmospheric parameters necessary for an 
abundance analysis of G and K giant stars by 
spectral synthesis or by equivalent widths and 
a fine analysis, both depending on good model 
atmospheres, has been described in several pa- 
pers (Luck, 1979; Luck and Bond, 1980). After 
a microturbulent velocity has been chosen from 
a curve of growth, the effective temperature is 
fixed at the value so that the Fe abundance 
derived from both Fe I and Fe II lines shows 
no dependence on the excitation potential of the 
lower level of the atom. The estimated error in 
the temperature determined by this means is 
about ±0.05 in 0 = 5040/7’ for early M giants, 
or about 150 K (Luck and Bond, 1980). No Fe 
II absorption lines appear in the spectra of M 
giants, however, and therefore the temperature 
determinations must be based solely on lines of 
Fe I. Observational evidence for departures 
from LTE in the K0 giant star, Pollux, has been 
presented by Ruland et al. (1980), and the ex- 
tra precautions against using lines from low- 
lying levels of Fe I is described by Luck (1982). 

A further instructive example is the detailed 
discussion of the derivation and uncertainty in 
effective temperature, surface gravity, abun- 
dance, and turbulent velocity for giant stars in 
globular clusters (Pilachowski et al., 1983), 
some of which are only slightly hotter than the 
red-giant stars of interest in this volume. The 
task is, of course, considerably easier in the 
globular cluster stars, for which a common 
mass (0.8 M© ) can be assumed with less resul- 


tant error in surface gravity (A log g = 0.25) 
than the value deduced from ionization 
equilibrium (A log g = 0.50). Departures from 
LTE appear to be the probable cause of the dif- 
ference in Fe abundance deduced from Fe I and 
Fe II. If so, the effect is an overionization of 
Fe I, the direction predicted theoretically 
(Auman and Woodrow, 1975). Observational 
evidence for overionization of Ca in M 
supergiants has been presented by Ramsey 
(1981). 

Effective temperatures for M giants derived 
from lunar-occultation angular diameters 
(Ridgway et al., 1980b) or from IR photometry 
(Tsuji, 1981a) and for M supergiants derived 
from a variety of sources have been plotted 
against the unreddened V-K color, and that 
relation has been used to infer the effective tem- 
perature for individual stars for which the (un- 
reddened) V-K colors can be obtained (Luck 
and Lambert, 1982). 

Thermal Structure 

Even when the effective temperature of a 
star is known, the thermal structure — T(t), 
T(P), or T(r ) — is a quantity of additional in- 
terest. In the Sun, this relation is easily deter- 
mined by limb-darkening measurements, but 
such measurements are essentially impossible 
in stars, except perhaps in special cases (cf. the 
discussion of limb darkening for Aldebaran by 
Ridgway et al., 1982). Possible methods for 
probing the structure of a stellar atmosphere 
are few: (1) a flux scan at a fixed wavelength 
across the solar or stellar disk, (2) a flux scan 
in wavelength from line center to wing across 
a strong spectral line, (3) observations at line 
center of different lines of a multiplet, (4) a flux 
scan in wavelength from the visible part of the 
spectrum toward shorter or longer wavelengths, 
and (5) the variation in angular diameter with 
wavelength. The fourth method probes increas- 
ingly shallower layers in the atmosphere be- 
cause the opacity increases in both directions 
due to atomic lines and Rayleigh scattering at 
shorter wavelengths and increasing free-free ab- 
sorption at longer wavelengths. 
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As far as the author is aware, no attempt 
has yet been made to deduce empirically, us- 
ing these or any other methods, the thermal 
structure of even a single red-giant star. All of 
our knowledge comes from the opposing direc- 
tion-comparing line profiles predicted from 
models with those observed — in which the fits 
can be said to be good (for example, see Quer- 
ci and Querci, 1975a; Olson and Richer, 1979; 
Johnson et al., 1982). However, these fits are 
not good enough, nor is the physics of line for- 
mation sufficiently well applied, to be used for 
diagnostics. We can perhaps say that the 
models are not badly off in thermal structure; 
as we have seen, however, even crude models 
of atmospheric structure and line formation 
yield profiles in good agreement with obser- 
vations. 

Ramsey and Johnson (1975) attempted to 
find sufficiently rich and strong multiplets to 
apply method (3) and invert their line-center 
fluxes to find a run of excitation temperature 
with depth (Ramsey, 1977). Unfortunately, the 
presence of both rotation and (large) macro- 
turbulence in real red-giant stars places severe 
restraints on the utility of this method (Desi- 
kachary and Gray, 1978). 

A study of the strong infrared lines of the 
first overtone bands of CO in Arcturus (K2 III) 
reveals an apparently irreducible discrepancy 
between models deduced from CO and those 
deduced from Ca II K, forcing a realization 
that inhomogeneous atmospheres are upon us 
(Ayres and Linsky, 1975; Heasley et al., 1978). 
An investigation of these lines in red-giant stars 
appears to hold considerable promise for de- 
ducing thermal structures. 

Of particular interest is the value of the tem- 
perature minimum. In the Sun, one sees radia- 
tion from the region of the temperature mini- 
mum in the ultraviolet (1600 A) and infrared 
(130 fim) wavelength regions (cf. Vernazza et 
al., 1981; Avrett, 1984). We have almost no 
knowledge of this quantity for any red-giant 
star except Betelgeuse (!T min = 2700 K), and 
even this value is uncertain because of the fail- 
ure of any homogeneous chromosphere to 
match both Ca II and Mg II line profiles (Basri 


et al., 1981). Although Johnson and O’Brien 
(1983) speculate that the spectral region near 
2850 A is formed at the temperature minimum 
in three N-type carbon stars, no confirmation 
by extensive line identification nor atmospheric 
modeling has yet been published. We anxiously 
await the first identification of the wavelength 
region corresponding to the temperature mini- 
mum in any red-giant star . 

Chemical Composition 

Much of the interest in red-giant stars arises 
from the unusual composition of their surface 
layers. Changes in the composition during the 
star’s lifetime are vital clues to the processes of 
nucleosynthesis and mixing which have taken 
place. These processes have been repeatedly dis- 
cussed from the point of view of the interior 
of the star (cf. Iben and Renzini, 1983, 1984). 
Our interest here is in inferring the abundances 
from stellar lines. In warmer stars, there is usu- 
ally no problem in finding sufficient atomic and 
molecular lines to identify and find the abun- 
dance of numerous elements, including even 
those rare species. In cool stars, the problem 
is much worse, and one must be careful to select 
those lines which give a valuable diagnosis of 
abundance. The lines and methods are carefully 
described in several papers (Wallerstein, 1973; 
Lambert, 1985; Gustafsson, 1985) to which 
readers are referred. 

Some of the difficulties involved in abun- 
dance analysis of the cool giants are described 
by Luck and Lambert (1982) in their attempt 
to extract the Li abundance in 3 1 M giants and 
supergiants. Even after careful calculation, in- 
cluding spectral synthesis of TiO lines over a 
1000 A spectral band, the level of the contin- 
uum near the Li I line remains the largest uncer- 
tainty in their analysis. 

The well-known difficulties in diagnosing 
the chemical composition of carbon stars are 
dramatically illustrated by V460 Cyg (earlier 
called DS Peg). Repeated analyses of high- 
dispersion spectra by empirical pseudocurves of 
growth (curves of growth using central depths 
rather than equivalent widths) and the isointen- 
sity method led to values of 12 C/ 13 C = 100 
(Fujita and Tsuji, 1977). Parallel analyses with 
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similarly high-resolution spectra using the line- 
depth method and empirical pseudocurves of 
growth by another group yielded 12 C/ 13 C = 8 
(Climenhaga et al., 1977). Later work aimed at 
reducing the discrepancy gave values of 32 from 
CN and CO (Dominy et al., 1978), of 25 from 
spectral synthesis of CN from the same spec- 
tra used by Climenhaga et al. (Johnson et ah, 
1982), and of 22 to 29 from the same method 
as Climenhaga et al. but with different spectra 
(Little-Marenin and Little, 1984). Recent work 
based on synthetic spectra and new model at- 
mospheres which include polyatomic opacities 
will provide more reliable values of CNO abun- 
dances and isotopic ratios in N-type carbon 
stars (Gustafsson 1985: private communica- 
tion). 

Surface Gravity 

Surface gravity not only is important in its 
own right but also is useful as a constraint on 
the mass, especially since stars of all masses 
tend to mingle in that part of the HR diagram 
where red giants are found. The emergent flux 
from a plane-parallel model is very weakly de- 
pendent on the gravity, and excellent observa- 
tions for stars for which other parameters (such 
as chemical composition and turbulent velocity) 
were well known would be necessary to give any 
leverage on the gravity. That the gravity of such 
a well-observed cool giant star as Arcturus 
should be sufficiently controversial to justify 
an international conference to investigate the 
matter testifies eloquently of our ignorance (cf. 
Trimble and Bell, 1981). The uncertainty in 
gravity increases as the temperature decreases. 

Gravities cannot be directly measured, but 
must be inferred from some other quantity, 
such as the pressure. A common method is to 
adopt for the gravity that value for which the 
abundance of a common element, such as Fe 
or Ti, is the same when two different stages of 
ionization are employed. Finding sufficient 
unblended lines from two such stages becomes 
increasingly difficult in the coolest stars. 

An illustration of the methods and problems 
in deducing the surface gravity comes from 
slightly hotter stars, in which a comprehensive 
scheme for deducing the effective temperature. 


surface gravity, and microturbulent velocity for 
G and K giant and supergiant stars has been 
outlined by Luck (1979). All of these parame- 
ters must be inferred from Fe lines because only 
that element is represented by a sufficient num- 
ber of suitable lines. The inferred abundance 
of Fe depends naturally on the value of all three 
of these parameters. One proceeds to choose 
that value of the microturbulent velocity for 
which there is no dependence on abundance for 
moderately strong lines (IF=200 mA). That ef- 
fective temperature is then chosen for which 
there is no dependence of the abundance on the 
excitation energy of the lower level of the line. 
That value of gravity is then chosen for which 
the abundance of Fe as deduced from Fe I and 
Fe II lines was the same. One then iterates as 
necessary, but this is not usually necessary. Sur- 
face gravities in K giant stars deduced by this 
scheme have estimated uncertainties of a fac- 
tor of 2. 

Unfortunately, because Fe II lines are weak 
or absent in the spectra of M giants and super- 
giants (except for the Fe II emission lines in the 
ultraviolet, which are probably circumstellar 
(Boesgaard and Boesgaard, 1976)), even the 
foregoing scheme is inapplicable, and to the au- 
thor’s knowledge, no application of the method 
of ionization equilibrium to red-giant stars has 
been made. In fact, some evidence cautions 
against it (Ramsey, 1981). The following pro- 
cedure has therefore been adopted by some 
workers (Luck and Lambert, 1982; Luck, 1982; 
and references therein). The luminosity of the 
star, generally estimated from the Wilson- 
Bappu effect and a bolometric correction, is 
combined with the effective temperature and 
compared with theoretical evolutionary tracks 
on an HR diagram, from which a mass is in- 
ferred. From the mass, effective temperature, 
and luminosity, the surface gravity is deduced. 
For M stars, the uncertainty in log g is 
estimated to be 0.5 corresponding to a factor 
of 3 in surface gravity or mass. 

An estimate of surface gravity can also be 
obtained from the pressure broadening of the 
wings of spectral lines — the traditional method 
in hotter stars. If the abundance is not accu- 
rately known (a common problem in red-giant 
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stars), that quantity must be determined from 
nearby weak lines formed in the same region 
of the atmosphere as the wings of the strong 
line to be used for gravity determination 
(Blackwell and Willis, 1977). The method, 
equivalent to finding lines on the damping por- 
tions of the curve of growth and calculating the 
damping constant, and hence the gravity from 
these, has not yet been applied to red-giant 
stars, but it should be useful if spectral lines 
meeting the necessary criteria can be found. 

A new method, specifically aimed at stars 
with molecular features, has recently been sug- 
gested and applied to Arcturus (Bell et al., 
1985). This method attempts to deduce the 
gravity from the pressure as the latter is mani- 
fested through molecular concentrations — in 
this case, MgH. Again, although this is not a 
cool star in the sense of the definition of this 
volume, the method has obvious application. 

Analysis of the Ca II K line core and wings 
holds potential for diagnosing the thermal 
structure, the turbulent velocity, and perhaps 
the pressure of cool outer atmospheres, but it 
has not yet been applied to red-giant stars. The 
method of formation of the line cores is still 
a matter of some uncertainty, with several con- 
tending theories in the field (cf. Linsky, 1980): 
(1) the width of the emission core (IF o ), mea- 
sured by Wilson and Bappu (1957) and cor- 
related with the absolute visual magnitude of 
the star, is formed in the Doppler core of the 
line and is therefore most sensitive to the turb- 
ulent velocities in the atmosphere; (2) the emis- 
sion core is formed in the damping part of the 
profile and hence is most dependent on the col- 
umn mass above the temperature minimum (cf . 
Ayres, 1979); (3) the line is formed in an iso- 
thermal layer (no chromosphere necessary) by 
the fact of a sharp drop in the source function 
as it changes to coherent scattering just outside 
the Doppler core (Basri, 1980); (4) the line is 
formed by a temperature rise due to an in- 
stability caused by infrared CO emission 
(Ayres, 1981; Kneer, 1983). Clearly the diag- 
nostic value of the line for various parameters 
depends on the mechanism of line formation, 
and application must await further elucidation 


of this mechanism, which may not be the same 
in every star. 

Is this relevant for red-giant stars? Can we 
expect the same processes to hold here? Several 
authors have derived empirical relations be- 
tween the width (FWHM) for the Ca K emis- 
sion and other characteristics of the star (e.g., 
one common formulation, based on 55 stars, 
is log W Q = -0.22 log g + 1.65 log T ff + 
0.10(Fe/H) - 3.69. Other authors find quite 
similar expressions (Ayres et al., 1975; Engvold 
and Rygh, 1978; Cram et al., 1979), including 
some based on theoretical calculations. Sikorski 
(1982) has generalized somewhat the relation 
found by Cram et al., and Glebocki and Stawi- 
kowski (1980) have extended it to include more 
fully the dependence of the line width on the 
core emission intensity. A similar relation holds 
for the Mg II h and k lines (Weiler and Oeg- 
erle, 1979). Unfortunately, the dependence on 
surface gravity is very weak and is of little value 
in stars in which the line itself is weak, as it 
tends to be in M giants. No application of these 
possibilities have yet been made in our stars ex- 
cept for Betelgeuse. The lines were measured 
in a few R stars and were used to find the 
luminosities (Richer, 1975) and thus confirm 
that the R stars were giants, whereas the hy- 
drogen-deficient R stars were supergiants. In 
fact, incredible as it may seem, no measure- 
ment, not even detection, of Ca K line cores 
has yet been made in cool carbon stars. 

Turbulence 

It is remarkable that, in a conference as re- 
cent and comprehensive as the IAU Collo- 
quium 5 1 on Stellar Turbulence (Gray and Lin- 
sky, 1980), hardly a word can be found on 
turbulence in red-giant stars. Understandably, 
theories must be tested under the simplest cir- 
cumstances and under conditions closer to 
those on the Sun, which is the standard astro- 
physical laboratory, but the paucity indicates 
our present poor knowledge of these processes 
in cool giants and supergiants. 

Methods for deducing turbulence are the 
classic ones, in which an extra broadening 
parameter beyond the thermal broadening is 
added to the line-absorption coefficient to 
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match the stellar line profiles. These methods 
are described under surface gravity because the 
two are inferred from the same data. Among 
early K stars, the macroturbulent velocities ap- 
pear to remain near 3 km s' 1 for stars with a 
range of M bol from +4 to -1, whereas for 
brighter stars, there is a sharp increase to about 
6 km s _1 but with considerable scatter (Smith 
and Dominy, 1979). Values are in the range 2.5 
± 0.5 km s' 1 for K and M giants (Luck and 
Lambert, 1982) to 4.0 ± 1.0 km s' 1 for M su- 
pergiants (Luck, 1982). Although values for in- 
dividual stars, especially some supergiants, are 
somewhat larger than these mean values, they 
are generally considerably lower than those 
given in the tabulation by Glebocki (1973), at 
least partly due to the lower value of turbulent 
velocity now accepted in the Sun. 

In a study of the lithium abundance in 19 
G and K lb stars, microturbulent velocities of 
1.9 to 2.6 km s' 1 and macroturbulent velocities 
of 4.0 to 9.0 km s' 1 were determined (Luck, 
1977a, 1977b, 1978). 

Very few reliable determinations of turbu- 
lent velocities in S or C stars have been carried 
out, but we note that, in a study of two N-type 
carbon stars by spectral synthesis and model at- 
mospheres, a microturbulent velocity of 4 to 5 
km s _1 and a weakly determined macroturbu- 
lent velocity of 5 km s _1 were required to 
match the broad lines of the CN 2-0 system in 
two N-type carbon stars (Johnson et al., 1982). 

With values of turbulent velocities are now 
available for several red-giant stars, these have 
been derived from a variety of methods and 
data and are not of the same precision. Tur- 
bulent velocities from Miras vary with phase, 
wavelength, and excitation potential because of 
their formation in different regions of an ex- 
tended and mass-losing atmosphere. Many of 
the velocities in the literature are in fact super- 
sonic; the sound speed in a red giant of T ff = 
3000 K is about 5 km/s, and it varies little for 
the conditions described here. Certainly no 
claim should be advanced for supersonic tur- 
bulence on the basis of any values inferred to 
the present, because deficiencies of various 
kinds in models and line formation can 


manifest themselves as microturbulence. There 
is also a downward trend of the values as spec- 
tra of higher resolution and better models have 
been employed. 

The turbulent pressure enters the equation 
of hydrostatic equilibrium through its contribu- 
tion to the total pressure. In a sense, it can be 
thought of as a levitating force slightly reduc- 
ing gravity. As far as the author is aware, no 
attempts have been made to deduce the effects 
of turbulent pressure in cool giants. Indeed, 
these depend on the same techniques as those 
for inferring the gravity and atmospheric ex- 
tension, and any information on turbulent 
pressure must await more precise measurement 
of these other quantities. 

Sphericity 

Although sphericity may be important in 
many of the coolest red-giant stars (especially 
in supergiants and Miras), insufficient knowl- 
edge of the influence of other factors makes it 
difficult to identify any observational features 
that are dependent primarily on extension. 
Lines of TiO in the red have been suggested as 
especially sensitive to atmospheric extension 
(Scholz and Wehrse, 1982). This area of re- 
search is just beginning. 

Chromospheres, Nonradiative Heating, 

Departures from LTE, and Mass Flow 

It was envisioned, when this volume was be- 
ing planned, that we would now proceed to a 
discussion of those spectral features inexplica- 
ble by thermal models, which, therefore, were 
diagnostic for the nonthermal and nonstation- 
ary processes that form the main thrust of this 
series of books (cf. Thomas, 1983). As has be- 
come painfully clear from at least this chapter, 
however, even common observations are still 
beyond the explicatory powers of present ther- 
mal atmospheric models. 

Quasi-thermal models, including chromo- 
spheres and departures from LTE, are analyzed 
in a succeeding chapter by de la Reza. As in K 
giant stars, the best diagnostic features for 
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stellar chromospheres are emission cores in the 
Ca II H and K lines, Fe II emission lines in the 
ultraviolet (Bidelman and Pyper, 1963; Boes- 
gaard and Boesgaard, 1976), emission or excep- 
tionally weak or strong absorption in Ha, and 
emission in Mg II h and k lines (cf. Linsky, 
1980, 1982; Baliunas, 1984; Brown, 1984). 
Lines from coronal and transition regions are 
generally not seen in red-giant stars because 
these stars all lie on the “mass-loss” side of the 
line which divides stars with coronae from stars 
with mass loss (Linsky and Haisch, 1979; Ayres 
et al., 1981). The prime diagnostic for mass loss 
is the violet-shifted central absorption compo- 
nent in the Mg II h and k lines (Stencel and 
Mullen, 1980). This latter requires a high-res- 
olution spectrum, and most types of red-giant 
stars are too faint for IUE. 

Spectra of M giant and supergiant stars al- 
most universally show all the foregoing indi- 
cators. However, except for the well-studied 
star Betelgeuse (Goldberg, 1984, and this vol- 
ume), models for the upper photosphere, the 
region of the temperature minimum, and the 
chromosphere are generally yet in the future (de 
la Reza, this volume). 

Much less is known about S stars because 
they have not been so intensively studied. An 
early study found variable Mg II h and k emis- 
sion in the S-Mira, Chi Cyg (Cassetella et al., 
1980), and emission in this line (h and k are 
blended at the low resolution needed to detect 
the ultraviolet emission in these cool objects) 
has been reported in several non-Mira S stars 
but with a wide range in strength (Johnson et 
al., 1984). No high-resolution ultraviolet spec- 
trum has been published. 

Certain N-type carbon stars have been 
known to show Fe II emission lines in the violet 
(Bidelman and Pyper, 1963), but no Ca H and 
K line data are available (Richer, 1975). Ultra- 
violet spectra taken by IUE show weak Mg II 
h and k lines in the few cool carbon stars that 
are bright enough to detect (Johnson and 
O’Brien, 1983). A single high-resolution spec- 
trum of TXPsc is now available and under 
study. The emission lines in HD 20234 are also 
variable (Querci and Querci, 1985). 


R stars show Ca II H and K emission cores 
(Richer, 1975), indicating a chromosphere, but 
low-resolution IUE spectra are remarkably sim- 
ilar to those of G and K giants. The early R 
stars (R0-R3) show, on low-resolution IUE 
spectra, Mg II and Mg I in absorption and no 
emission features whatever. Middle R stars (R5) 
show either no Mg II h and k line or faint ab- 
sorption features. Late R stars (R8) show weak 
Mg II emission features similar to those of N- 
type stars (Eaton et al., 1985). Only one high- 
resolution IUE spectrum has been obtained (the 
R0 HdC star, HD 182040), and this shows only 
the faintest hint of any emission at h or k. 

Chromospheric modeling for red-giant stars 
had hardly begun (cf. Avrett and Johnson, 
1984) and presents a perfect illustration of a 
field ripe for research (de la Reza, this volume). 
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